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43 ! Abstract 
Oh! 

We prove that Euclidean 4 -quantum field theory on four-dimensional 
Moyal space with harmonic propagation, taken at its self-duality point 
and in the infinite volume limit, is exactly solvable and non-trivial. 
The theory is formulated as a matrix model with partition function 
Z[J] = J dM exp (trace (JM - EM 2 - f M 4 )). The integral is over the 
space of Hermitean M x AAmatrices, E is an external matrix which for 
J\f — > oo has the same spectrum as the Laplace operator in 4 dimensions, 
A > is a scalar coupling constant and the matrix J is used to generate 
^ | correlation functions. 

Using Schwinger- Dyson techniques and the theory of singular integral 
CNj ■ equations of Carleman type we compute (for Af — > oo and after renor- 

malisation of E, A) the normalised partition function Z[J]/Z[0] exactly 
in terms of the solution of a non-linear integral equation. Existence of 
a solution is proved via the Schauder fixed point theorem. The deriva- 
tion of the non-linear integral equation relies on an assumption which we 
verified numerically for coupling constants < A < 647r. 
Correlation functions for — > oo have the topology of a 2-sphere with B 
punctures. We prove an algebraic recursion formula for all functions with 
a single puncture and list the solution for the 4-, 6- and 8-point functions 
at B — 1 explicitly, thereby showing that these functions are non-trivial. 
Bare and effective coupling constants differ only by a finite ratio, which 
means that the beta-function is non-perturbatively zero. The (1+1)- 
and (2+2)-point functions, and conjecturally any (Ni+ . . . +A ? B )-point 
function with iVj < 2, are given by a more complicated but still explicit 
formula. Any function with at least one iVj > 3 is again algebraic. 
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1 Introduction 



A rigorous construction of quantum field theories in four dimensions was not very suc- 
cessful so far. In this paper we show that for 4 -theory on four- dimensional Moyal space 
with harmonic propagation, taken at critical frequency and in the infinite volume limit, 
much more is true: The model is exactly solvable. 

We know that this is a toy model in so far as classical locality and Poincare symmetry 
are not realised and the Minkowskian continuation of that Euclidean model needs to be 
investigated. On the other hand, the model 

- carries an action of an infinite-dimensional symmetry group, 

- is invariant under a duality transformation between position and momentum space 
ESU2] , 



is almost scale invariant IDGMR07 



is known to have a realisation as a matrix model (with non-constant kinetic term), 
has perturbatively an infinite number of divergent but renormalisable Feynman 



graphs ||GW05b 



Each renormalised Feynman graph has subleading logarithmic terms which make the 
perturbatively renormalised correlation functions divergent at large energy. The model 
and its solution touch many aspects of quantum field theory which we recall in the next 
subsections. The reader in a hurry may jump to Sec. [O. 



1.1 Perturbative, axiomatic and algebraic quantum field theory 

Starting with the Lamb shift in the 1940s and culminating in the experimental tests of 
the Standard Model, perturbatively renormalised quantum field theory is an enormous 
phenomenological success. However, this success lacks a mathematical understanding. 
The perturbation series is at best an asymptotic expansion which cannot converge at 
physical coupling constants such as the electron charge e 2 In addition there are 

physical effects such as confinement which seem out of reach for perturbation theory. 
Therefore, the development of a mathematical foundation of quantum field theory that 
permits non-trivial examples is one of the most urgent challenges in mathematical physics. 

In the early 1950's, Garding and Wightman gave a rigorous mathematical founda- 
tion to relativistic quantum field theory by casting the unquestionable physical principles 
(locality, covariance, stability, unitarity) into a set of axioms. These ideas were published 
years later [|Wig56| , |WG64| , |SW64|1 . The difficulty to provide non-trivial examples to these 



axioms motivated the development of equivalent or more general frameworks such as Al- 
gebraic quantum field theory and Constructive/Euclidean quantum field theory. Algebraic 
quantum field theory shifts the focus from the field operators to the Haag-Kastler net of 
algebras assigned to open regions in space-time ||HK64|| . Fields merely provide coordinates 
on the algebra. Over the years this point of view turned out to be very fruitful |)Haa92 . 



One important advantage over the axiomatic setup is the natural possibility to describe 
quantum field theory on curved space-time [|BF09| . 
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1.2 Euclidean quantum field theory 



Of central importance for us is the Euclidean approach. Wightman functions admit 
an analytic continuation in time. At purely imaginary time they become the Schwinger 
functions [ fSch5!J ] of a Euclidean quantum field theory. Symanzik emphasised the powerful 
Euclidean-covariant functional integral representation ||Sym64|| , which yields a Feynman- 
Kac formula of the heat kernel ||Kac49|| . In this way the Schwinger functions become the 
moments of a statistical physics probability distribution. This tight connection between 
Euclidean quantum field theory and statistical physics led to a fruitful exchange of con- 
cepts and methods, most importantly that of the renormalisation group | WK74 . 

It is sometimes possible to rigorously prove the existence of a Euclidean quantum field 
theory or of a statistical physics model without knowing or using that this model derives 
from a true relativistic quantum field theory. This is, for instance, the case for the model 
constructed in this paper. Sufficient conditions on the Euclidean model which guarantee 
the Wightman axioms were first given by Nelson [|Ncl71||. These conditions based on 



Markov fields turned out to be too strong or inconvenient. Shortly later, Osterwalder and 
Schrader established a set of axioms ||OS73|1 by which the Euclidean quantum field theory 
is equivalent to a Wightman theory. The most decisive axiom is reflection positivity which 
yields existence of the Hilbert space and a positive energy Hamiltonian. The Euclidean 
approach together with the Osterwalder-Schrader axioms turned out to be the key to 
construct relativistic quantum field theories in dimension less than four. Two successful 
methods were developed: The correlation inequality method ||Gin70| , |GRS73|| relies on 



positivity and monotonicity of Schwinger functions and is suitable for bosonic theories. 
The expansion method [ |Riv91|| establishes bounds for the coefficients in the perturbation 
series so that for small positive coupling constants the series can be resummed. We also 
mention the cluster expansion ||GJS74|| which is a key to prove a mass gap in the spectrum 
of the Hamiltonian. 



1.3 Solvable models 



Under solvable models we understand models in quantum field theory or statistical physics 
where all correlation functions can be exactly evaluated in terms of "known" functions, 
at least in principle. With a few three-dimensional exceptions, solvable models were 
established only in one or two dimensions. 

The first field-theoretical example was the Thirring model ||l'hi58|| which describes a 
quartic self-interaction of a Dirac field in 1+1 dimensions. In the massless case, Johnson 
||Joh61|| found the exact expression for the 2- and 4-point functions. This work was 
extended by Hagen ||Hag67|| and Klaiber [|Kla67|| to the explicit solution of any correlation 
function. The massive case is more complicated; we mention the construction of the S- 
matrix by Korepin | |Kor79| using the Bethe ansatz | [Bet31| . Another famous model that 
is exactly solvable is the Schwinger model ||Sch62|| , or QED in 2 dimensions. 

Many 2-dimensional models of statistical physics are known to be solvable. The first 



one was the 2-dimensional square-lattice critical Ising model [Isi25], which after important 



work by Kramers and Wannier [KW41| was solved by Onsager ||Ons44|| using the transfer 
matrix method. The list of more involved exactly solvable models contains the 6-vertex 
model (or ice model) solved by Lieb [|Lie67|l , the 8- vertex model solved by Baxter |Bax71 



and the hard hexagon model also solved by Baxter ||Bax80|| . The quantum inverse scat 
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tering method | |Fad95| | and Yang-Baxter equations | JJ im9(J| ] are important tools for these 
achievements. 

The Ising model turned out to be the prototype of a discrete series of solvable two- 
dimensional field theories: the Minimal Models in conformal field theory [ |BPZ84| , |FQS84 



Minimal Models correspond to highest-weight representations of the Virasoro algebra 

at central charge c = 1 , 6 , ,s . Some of these models give non-trivial realisations 

of the Osterwalder-Schrader or Wightman axioms. A solvable model of different type 
in conformal quantum field theory is the Wess-Zumino-(Novikov-)Witten model ||WZ71 



Wit84|| whose solutions are realised by affine Kac-Moody algebras. 



1.4 Constructable models 

Constructing a model means to prove for a specific candidate Hamiltonian or Lagrangian 
the axioms of Wightman, Haag-Kastler or Osterwalder-Schrader. Thereby the Wightman 
or Schwinger functions, although not computed/solved explicitly, are shown to have the 
required properties. Constructing a model is less than solving it. However, most of the 
solved models are disguised free fields whereas the constructed models are true interacting 
field theories, unfortunately only in dimension < 4. For a historical review (that we used 
already in the section on Euclidean quantum field theory) we refer to | JafOO|| . 

The first successfully constructed model is the 4 - model in 2 dimensions, §\ for short, 
constructed in a series of articles by Glimm-Jaffe | GJ68-72" |. Whereas this was directly 
achieved in the Wightman or Haag-Kastler setup, the Euclidean approach was decisive for 
the generalisation to the P{<p)2 model ||GJS74j , |GRS73|| , i.e. a two-dimensional scalar field 
with polynomial interaction. The construction of <p\ [ GJ73 . FQ76 | turned out to be much 



harder. There is a stability problem which was overcome by renormalisation group ideas. 
Another problem not present in two dimensions is the inequivalence of representations of 
the canonical commutation relations for the interacting and the free field. 

Perturbatively, these models 02, 03 and P{4>)2 are all super-renormalisable. It is 
therefore not surprising that a construction of the 04-model (which is perturbatively just 
renormalisable) faces far more problems. In 4 + e dimensions, the 4 -model does not exist. 
As shown by Aizenman and Frohlich |Aiz81| , [Fro82|| , the model is trivial , i.e. it only exists 
if the coupling constant vanishes. Related to triviality is the appearance of Landau poles 
| LAK54j| , first discovered for quantum electrodynamics (QED 4 ): The /3-function which 
describes the running of the coupling constant when changing the scale could develop a 
singularity at finite momentum cut-off. This means that the theory is only consistent 
below that cut-off, or has to be trivial if the cut-off is removed. Perturbatively, the /3- 
function of 04-theory develops a Landau pole, which together with the triviality of 
strongly suggests (although a proof is still missing) that 



b\ cannot be constructed. 



There is a variant of the Euclidean A^j-model, the wrong-sign model A < 0, which 



can be constructed if restricted to planar graphs ||tHo82| , |Riv83|| . Because of the wrong 



sign and the resulting instability, this Euclidean model does not give rise to a model of 
relativistic quantum field theory. We mention this model because it shares some aspects 
with the model we study and solve in this paper. Our model also restricts to the planar 
sector, but as result of the thermodynamic limit and not by hand. Our model is stable 
but violates locality in the traditional sense. It is unknown so far whether or not our 
model satisfies a non-local variant of the Osterwalder-Schrader axioms. 
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Recall that QED4 is ruled out for the same Landau pole problem [|LAK54 |. It is there- 
fore important that non-abelian Yang-Mills theory is asymptotically free ||GW73| , Pol73 
and as such does not have a Landau pole (in perturbation theory). This means that 
Yang-Mills theory is a candidate for a constructable quantum field theory in four dimen- 
sions. Unfortunately, Yang-Mills theory is too difficult so that its construction is an open 
problem [ jj W00| |. A toy model of asymptotic freedom is the Gross- Neveu model [ GJN 74 1 , 
a generalisation of the Thirring model to A^-component fermions. In two dimensions this 
model (GN 2 for short) is perturbatively renormalisable but not super-renormalisable. The 
construction of (GN2) was achieved in ||GK85| , [FMRS85| and constitutes the first example 
of a constructable just renormalisable quantum field theory. 



1.5 Noncommutative geometry 

Our present fundamental physics rests on two pillars: Quantum field theory and General 
relativity. One of the main questions in this area concerns the matching of these two 
concepts. In perturbative quantum field theory, the standard model which correctly 
describes all experimentally observed particles is renormalisable, whereas general relativity 
which correctly describes all observed gravity effects is not. The renormalisation group 
tells us that non-renormalisable interactions are scaled away. To put it differently, the 
existence of gravity means that space-time cannot be viewed as a differentiable manifold 
over all length scales. From the numerical value of the gravitational coupling constant the 
renormalisation group estimates the fundamental scale at which the manifold structure 
breaks down: it is the Planck scale 10 -35 m. The most active approaches of quantum 
gravity all agree on the assumption that space-time is fundamentally different at the 
Planck scale; there is only disagreement about the structure which replaces it. 

The approach to quantum gravity which is relevant for this paper is noncommutative 
geometry | |Gon94| | . Noncommutative geometry has its roots in the mathematical descrip- 
tion of quantum mechanics. It has been vastly developed over the years and achieved 
spectacular success in mathematics. Its relevance to physics comes from the fact that both 
Yang-Mills theory (as a classical field theory) and general relativity are unquestionably 
of geometrical origin. Noncommutative geometry achieved, on the level of classical field 
theories, a true unification of Yang-Mills theory, and even of the whole standard model, 
with general relativity ||Con96| , |CC96|| : General relativity put on a noncommutative space 
which is the product of a manifold with a discrete space is nothing but Yang-Mills theory. 
The astonishing picture which results is that the breakdown of the manifold structure of 
space-time already occurs at the much larger length scale given by the Compton length 
of the W and Z bosons of the standard model. This breakdown is mild; what at larger 
distances is a point of space-time becomes a couple of points at scales shorter than the 
standard model scale. But once accustomed to this idea it is natural to assume that 
passing to shorter and shorter distances there could be a cascade of phase transitions in 
which the manifold structure of space-time fades out more and more. 

Quantum field theory as defined by the Wightman, Haag-Kastler or Osterwalder- 
Schrader axioms relies on the manifold structure of space-time. In giving up the manifold 
one has to adapt the axioms. A reasonable replacement is not yet available. One of 
the reasons is that whereas noncommutative geometry has a clear concept of compact 
manifolds with Euclidean signature [|Con08||, it is unknown how to describe non-compact 
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manifolds with Lorentzian signature. Therefore, the focus has been on Euclidean quantum 
field theories on noncommutative manifolds. Because of their functional integral realisa- 
tion, such models are easy to define: It suffices to specify a parametrisation of the fields 
and an action functional for them. Scalar fields, for instance, can classically be viewed 
as sections of a vector bundle over the manifold. In noncommutative geometry, such sec- 
tions form a projective module over an algebra, with the algebra itself being the simplest 
example. The traditional pointlike interaction of scalar fields is then translated into the 
product in a noncommutative algebra. 

The simplest examples are given by deformations. This means that one takes the same 
continuous or smooth functions on the manifold as parametrisation but equipped with a 
deformed product *. Rieffel showed [|Rie93| that if the manifold carries an action of M. d by 



translations, as it is the case for Euclidean space, it is possible to turn this group action 
into a deformed associative product. Doplicher-Fredenhagen- Roberts studied quantum 
space-time arising from space-time uncertainty relations at the Planck scale ||DFR95|| , 
Filk derived the Feynman rules for deformed Euclidean space |Kil96|| , and in 1999 many 
authors showed perturbative one-loop renormalisability of a number of models. Shortly 
later Minwalla, van Raamsdonk and Seiberg discovered a severe problem at higher loop 
order, the so-called ultraviolet /infrared mixing ||MVS00|| . 

A possible way to cure this problem for the 0|-model has been found by us in previous 
It leads to an action functional which has 4 relevant /marginal operators 



work 

Z, /i 2 , A, Q (instead of 3): 



S 



dXy — 



H —<P*< 



x 



Here is the deformation matrix which defines the ^-product. We have been able to 
show that the resulting model is renormalisable up to all orders in perturbation theory. 
In addition a new fixed point appears at O = 1. At one-loop order the theory flows 
into this fixed point, leaving the ratio ^ constant ||GW04| , |GW05cj| : In contrast to the 



usual 04-theory, the running coupling constant is one-loop bounded. The group around 
Rivasseau then showed that for the critical case Q = 1 the /3-function vanishes up to three 
loop order |DR(J7|| , and finally succeeded in proving (3 = up to all orders in perturbation 
theory PGMR07[| . In this way, the Landau pole problem is tamed. 

Vanishing of the /3-function is often a consequence of integrability. It was therefore con- 
jectured that the model ( |1.1| ) at critical frequency Q = 1 has a chance to be constructable. 
The construction was searched in two directions. There is the loop vertex expansion de- 
veloped by Rivasseau [|Riv07b | which combines the Hubbard-Stratonovich transform with 
the Brydges-Kennedy-Abdesselam- Rivasseau forest formula |BK87| , |AR94|| . This method 
already proved useful in traditional constructive quantum field theory [ |MR08| . |RW 1 2|| and 



has been applied in ||Wanl2|| to the two-dimensional version of (|1 . 1|) . The second approach 
was proposed by us in | jGW09| |. By extending the Ward identities and Schwinger-Dyson 
equations used by Disertori-Gurau-Magnen- Rivasseau in IPGMROTj] , we derived a self- 
consistent non-linear integral equation for the renormalised 2-point function of the model 



( p..l|) . This is a non-perturbative equation, which we solved perturbatively up to third 
order in the coupling constant. 
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1.6 Matrix models 



The large- M limit of SU(N) Yang-Mills theory led to matrix models with dominant 
planar diagrams [|tHo74|| . Since that observation, matrix models play an important role in 



mathematical physics: Triangulations of two-dimensional manifolds as simple models of 
quantum gravity coupled to matter led to many studies of matrix models ||DGZ95|| . The 
IKKT [[IKKT97I1 and BFSS [|BFSS97j1 models are important as reductions of string theory 
models. Kazakov et al [[BK90|| observed that in the double scaling limit of matrix models 
a new phase transition occurs, which is of great interest. 

The partition function of a one-matrix model as toy model for 2D gravity is given 
by Z = J dM exp(— a n tr(M™)), where the integral is over (A/" x A/j-Hermitean 
matrices and the a n are scalar coefficients which may depend on N '. In the double scaling 
limit a n = Mt n and Af — > oo, the partition function becomes a series in (t n ) which 
can be expressed in terms of the r-function for the Korteweg-de Vries (KdV) hierarchy 
PGZ95 



sec 



for a review). There is another approach to topological gravity in which 
the partition function is a series in (t n ) with coefficients given by intersection numbers 
of complex curves. Witten conjectured |[Wit91|| that the partition functions of the two 



approaches coincide. This conjecture was proved by Kontsevich ||Kon82|| who achieved the 
computation of the intersection numbers in terms of weighted sums over ribbon graphs 
(fat Feynman graphs), which he proved to be generated from the Airy function matrix 
model (Kontsevich model) 



/ dM exp ( - \ti(EM 2 ) + 

Z[E] = I 

/ dM exp ( - ±tr(£M 2 )) 



itr(M 3 )) 



1.2) 



where i 2 = —1 and E is a positive Hermitean matrix related to the series (t n ) by t n = 
{2n — l)!!tr(£ ,_ ( 2n_1 ^). The large-A/" limit of ( |1.2| ) gives the KdV evolution equation, thus 
proving Witten's conjecture. 

The Rieffel deformation of Euclidean space can be described by a matrix product 
GV88|| . Langmann, Szabo and Zarembo showed ||LSZ04|1 that the model for a complex 
scalar field $ with interaction ($ * $) 2 in an external magnetic field B = @ _1 which is 



dual ||LS02|| to the deformation matrix gives rise to a matrix model with partition function 
J dMdM ] exp ( - wi. EMM - \M MM M)). They showed that in the large- M limit 
the model is exactly solvable but trivial. In joint work with Steinacker ||GS06| , pS08 



one 

of us used the relation of the Langmann-Szabo self-dual noncommutative 0* 3 -model to 
the Kontsevich model to perform a non-perturbative renormalisation and to extract the 
running of the asymptotically free 0* 3 -model in 6 dimensions. 

In this paper we relate the Langmann-Szabo self-dual 0* 4 -model (which in contrast to 
<f)* 3 is stable) in 4 dimensions to a matrix model of partition function 



Z[E,J,\] 



J dM exp ( - tr(£M 2 ) + tr (JM) - f tr(M 4 )) 
dM exp ( - ti(EM 2 ) - \ tr(M 4 )) 



1.3) 
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The matrix E is fixed by the identification with the noncommutative field theory, and the 
external source J is used to generate the correlation functions. We identify a scaling limit 
M — > oo in which the model is exactly solvable. We are able to compute all correlation 
functions in terms of a single function G of one variable. The correlation functions show an 
intriguing combinatorial structure. In contrast to the Langmann-Szabo-Zarembo model, 
our model is non-trivial. It is tempting to speculate that there might be relations to 
mathematical structures which are the analogues of KdV hierarchy and intersection theory 
of complex curves to which the Kontsevich model is related. 



1.7 Outline of the paper 



This paper is a far-reaching generalisation of our previous work 



We first define 



in Sec. ^| a general framework of field-theoretical matrix models for a compact operator 
on a Hilbert space. Whereas the interaction is, as usual, the trace of a polynomial 
in 0, the kinetic term is tr(E<f) 2 ) for a selfadjoint unbounded operator E with compact 
resolvent. The action of the group of unitaries on Hilbert space induces an infinite number 
of Ward identities. In Theorem ^]3| we turn this Ward identity into a formula for the 
second derivative of the partition function. This is the most crucial step which goes far 



beyond previous glimpses | DGMR07 . |GW09 l of that formula. It is then straightforward 
to algebraically derive in Sec. ^ the Schwinger-Dyson equation of any correlation function 
(i.e. including non-planar functions and with several boundary components) of the matrix 
model, and not only the planar 2- and 4-point functions as in [|GW09 1 . For the planar 
regular 2-point function this results in a non-linear equation for that function alone. In 
sharp contrast with other models there exists an ordering of all other correlation functions 
such that the equation for a given function only depends on functions of smaller order. 
Except for the planar regular 2-point function, these equations are only linear in the top 
degree function and therefore "easy" to solve. This means that solving the non-linear 
equation for the planar regular 2-point function (which is hard) amounts to solve the 
entire model. 

So far this applies to any field-theoretical matrix model with cut-off. One example 
is the model (|1.1|) which we recall in Sec. f|. By results of Gracia-Bondia and Varilly 
||GV88|| , the model has a matrix realisation. We perform the mass renormalisation and 
prepare the wave function renormalisation. As shown by us and together with Gayral 
in |GW12| , |GW11| , the model has, as a noncommutative geometry, a finite volume of 

diameter \J\, where 9 is the deformation parameter which determines the ^-product. 

Since fi = 1 is fixed, going to the thermodynamic (infinite volume) limit amounts to 
the limit 9 — > oo. The limit (Q = 0, 9 — > oo) was already emphasised by Minwalla, 



Van Raamsdonk and Seiberg in ||1V1VS0(J|| and then further studied by Becchi, Giusto and 
Imbimbo in ||BGI02| . PGI03|| in the renormalisation group framework. Becchi et al found 
for graphs in momentum space, and we confirm this for whole functions at Q = 1, that 
the non-planar sector is scaled away but there remains a non-trivial topology in form of 
spheres with several holes (they called them 'swiss cheese' in [BGI02]| ). 

Moreover, and this is the key to all further progress, the infinite volume limit 9 — > 
oo turns the Schwinger-Dyson equation for the regular (B = 1 hole/puncture) 2-point 
function into an integral equation (Sec. ^J). Performing the wavefunction renormalisation 
in order to prepare for the continuum limit A — > oo, this integral equation for the 2- 
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point function G ab at continuous "matrix indices" a, b e [0, A 2 ] is cubic in G a b and has 
a singular integral kernel. In our previous work ||GW09|| we have already arrived at this 
point. The focus there was put on a perturbative solution up to third order in the coupling 
constant which revealed certain number-theoretic properties. In this paper we arrange 
the equation in different manner. In this way we recognise that the equation for G a b is 
linear in the function with a, b 7^ and non-linear only in its boundary value G a0 . Even 
more, the resulting linear equation at given boundary G a0 is the well-studied Carleman 
type singular integral equation | |Car22| , [lri57| , [Mus65| |. Using the known solution of that 



particular Riemann-Hilbert problem, we solve exactly the equation for G a b in terms of 
G a o. But since the 2-point function is symmetric, and Gob = Gbo is solved in terms of 
G a o, this gives rise to a consistency equation G = TG 

Gbo = t^t exp f - -A- tdt f°°dp — - 2 j | . (1.4) 

Here, is the Hilbert transform U™\f(*]\ = - lim lim ( ! + [ )dp 

TT A^oo e^O \ J Q J a+e ) p-Q 

The big progress over ||GW09|| is that the non-linear operator T is, although still 



complicated, an improvement operator: We prove that for A > the operator T satisfies 
the assumptions of the Schauder fixed point theorem which therefore guarantees that 
( |1.4| ) has a solution {b \-t Gbo) £ Cq(M + ) which is monotonously decreasing, positive, and 
vanishing together with its derivative at infinity. 

To be precise, the consistency equation ( fL.4| ) relies on an assumption which we do not 
check. In an accompanying paper ||GW13|| we study (|1.4|) numerically and validate this 



assumption for coupling constants < A < A* = 647r. This is far into the strong-coupling 
regime. For A > A* we find that the required symmetry G a b = Gb a does not hold for the 
numerical solution of Ql.4j) . We trace this back to the non-trivial solution |Tri57|| of the 



homogeneous Carleman equation which we have put to zero in the derivation of ([l^ 
The additional term is related to a winding number [|Mus65|| for the boundary value of 



an analytic function. It is very well possible that at A* one enters another sheet of the 
logarithm which would provide a continuation of the solution to A > 6Att. Moreover, the 
non-trivial solution of the homogeneous Carleman equation also affects the uniqueness 
proof. All this is left for future investigation. 



In sec. we fulfil our promise given in |GW09] and in the first version of this paper 



by showing how the Schwinger-Dyson equations for all other correlation functions can 
be solved in terms of the "known function" Gbo- For Appoint functions which have the 
topology of a sphere with B = 1 puncture we derive an algebraic recursion formula in A" 
and give the explicit solution for A" e {4,6,8}. The solution for the 4-point function is 

remarkably simple, G bobl b 2 b 3 = ^ ^1%^^ ' where C 1 + ^ = "T^U" We 
prove that the Appoint function has a similar structure, namely a sum of fractions where 
the numerator is a non-crossing complete pairing of the indices bo, ... , b^-i put in cyclic 
order on a circle and the denominator a product of a tree connecting all even indices with 
a tree connecting all odd indices, subject to the condition that any of the trees intersects 
the pairings only in the vertices. The theory is non-trivial, but easy enough to make 
this explicit solution possible. At coinciding indices, such as for the effective coupling 
constant \ e ff= —Gqqqq, one encounters derivatives of the 2-point function. We determine 
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X e ff explicitly in terms of A and the fixed point solution Gbo and find that X e ff differs only 
by a finite factor from the bare coupling A. We thus establish at non-perturbative level 
the perturbative result of Disertori, Gurau, Magnen and Rivasseau |DGMR07| that the 
model has vanishing ^-function. 

In the appendix we derive the Schwinger- Dyson equations for all (Ni+N 2 ) -point func- 
tions which have the topology of a sphere with two punctures. The (l+l)-point function 
G a \ c and the (2+2)-point function G a b\ c d are more involved but still computable in terms 
of Gbo- Similar to the iV-point functions, all (N1+N2) -point functions with iVj > 3 are 
sums of difference quotients of polynomials in G a b and G a \ c (for Ni odd) or G a b\cd (for Ni 
even). Although we do no work out the details, it is clear that any correlation function 
with topology of a sphere with B > 2 punctures can be computed along the same lines. 
This means that solving the fixed point equation (p..4|) amounts to solve the whole theory. 
This solution is a sum of infinitely many renormalised Feynman graphs. We check in 
appendix |C| that the perturbative expansion (at next-to-lowest order) of our exact results 
coincides with the Feynman graph computation. 

We stress again that this is an exact solution of a Euclidean quantum field theory. But 
thanks to the explicit knowledge of all correlation functions it should be possible in not 
too distant future to decide whether or not the model satisfies reflection positivity. An 
affirmative answer is not impossible in view of the positivity and monotonicity properties 
of the 2-point function. We are convinced that in perturbation theory there is no chance 
to settle this question. We feel that getting thus far with this model was only possible 
because there is a deep mathematical structure behind which is still to be discovered. We 
hope that the novel methods we developed and the fascinating results we obtained might 
be relevant also for realistic quantum field theories in four dimensions. 



2 Ward identities in matrix models 

A classical scalar field on M. d is (e.g.) a continuous function <p G Co(M d ) vanishing at 00. 
The space of such functions forms a commutative C*-algebra and hence is realised as a 
bounded linear operator on Hilbert space. One natural translation of existence of the 
mass term J R ddx ^-0 2 (x) to noncommutative operator algebras is to require that 2 is 
a trace-class operator, which is the case if is a Hilbert-Schmidt compact operator on 
Hilbert space. 



2.1 Field-theoretical matrix models 

A matrix for us is a compact operator on Hilbert space H. For the Hilbert space H = C n 
this is an ordinary n x n- matrix, for H = £ 2 (N) an infinite matrix and for H = L 2 (M. n ) 
an integral kernel operator. Accordingly, we let / be a set of indices, which can be finite, 
countable or continuous. 

We let Ai be the algebra of matrices = (0 a fe)a,fee/ associated with /, i.e. the vector 
space of maps / x I 3 (a,b) !->■ O & G C equipped with the product (4>ip) a b = J2 c &i 4><u$do- 
For continuous / the space Ai is the space of integral kernels (p a b = 4>{ a i fy, and the above 
sum is a short-hand notation for the integral J2cei 4>acipcb = f I dfJ,(c) (f>(a,c)i/;(c,b) with 
respect to a measure \l on /. It might be convenient to also allow for weight factors in 
discrete sums, (0V)afe = J2 ce i Vc<l>aci>cb- 
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We recall that for continuous / a matrix (0 a &) defines a compact operator on H = 
L 2 (I, /i) (in fact a Hilbert-Schmidt operator) if (0 a &) G L 2 (J x /, /i x /i). The case of 
infinite discrete I is then a special case for the Dirac measure. 

For tp = (ip a b) in the trace-class ideal of compact operators we let tr(-0) := XLVVa be 
the trace (tr(^) := Ylia^a^aa m presence of a measure), and 0* b = 0& Q defines the adjoint 
matrix of 0. 

We consider Euclidean quantum field theory for a single real scalar field = 0* G Ai 
defined by an action functional 

S[0]=tr(P0 2 + P[0]). (2.1) 

Here, P[0] is a polynomial in with coefficients in K, and V[<f>] = tr(P[0]) would be a 
standard matrix model action. Being interested in a situation closer to field theory, we 
need an analogue of the kinetic term J Rd dx ^<j)(x)(—A<p)(x) involving the Laplacian. We 
therefore let E to be an unbounded selfadjoint operator on Hilbert space with compact 
resolvent. This means that the resolvent (E — i) _1 is a compact operator on H, i.e. a 
matrix, so that by little abuse of notation we can view E as a matrix, too. To make sense 
of tr(P0 2 ) we require to belong to an appropriate subspace of the compact operators 
such that P0 2 is trace-class. Note that the Kontsevich model ||Kon82|| with action (|1.2|) 
fits into this setting, with the resolvent of E being trace-class. 



The action (|2.1|) gives rise to the partition function 

Z[J] = J V[<t>] eMS[<t\ +tr(0J)) , (2.2) 

where T>[4>] = Yl a b<=i d<t>ab is the extension of the Lebesgue measure defined on finite-rank 
operators to the space of selfadjoint matrices <fi G Ai for which E<p 2 is trace-class. By 
J G Ai we denote the external source. The partition function gives rise to connected 
correlation functions by 

(Va lbl ■ ■ .<Pa nbn ) ■= « T dnW[ i\ , W[J] := logZ[J] . (2.3) 



Unless / is finite, the resulting index sums may diverge and require a renormalisation. 
2.2 Ward identity 

The action of the group U(H) of unitaries on Hilbert space induces an adjoint action of 
U G U(H) on bounded operators, Ai 3 <p ^ U<f)U* G Ai. We let <fi u := U<pU*. Since the 
adjoint action preserves the space of selfadjoint matrices and is a dummy integration 
variable in the partition function, one has 



Z[J] = J V[<j>) exp(-5[0] +tr(0J)) = / V\<f\ exp(— 5[0 ] + tr(0 J)) . (2.4) 

For finite matrices A = A*,B G M n (C), the Jacobian of the transformation A i-> BAB* 
is det^^ = (detPdetP*)™. In particular, det ^ff- = 1 for U G U(n). Therefore, 
the transformation formula for [0] yields 

= f V[4>] { exp(-S[0] + tr(0J)) - exj>(-S[U(f)U*) + tr(U<f>U*J))\ (2.5) 
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for any unitary U. Note that the integrand itself does not vanish because tr(£'0 2 ) and 
tr(<pJ) are not invariant under i— > <p u . Linearisation of (|2.5|) about U = id gives for the 
action S[</>] = tr(E(p 2 + V[<f>]) the system of Ward identities 

= J V[4>] (E<p 2 - <P 2 E -J<p + 0j) exp ( - S[4>] + tr(0J)) . (2.6) 

Note that Q2.6| ) is a matrix equation, which means that for an infinite index set / there is 
an infinite number of (scalar) Ward identities. By adjoint action with U we can always 
achieve that the external matrix E is diagonal, E mn = E m d~ mn . In this case the system 
takes the form 

0= V[(j>] ^ (i E P ~ E a)4>pn4>na ~ Jpn4>na + Jna^pn^ exp ( - S[(j>] + tr(0j)) . (2.7) 



Writing under the integral (f) a b = g§—, we finally obtain (this was first derived in 



| DGMR07| ): 



Proposition 2.1 The partition function Z[J] of a real scalar matrix model with action 
S[4>] = tr(E(p 2 + P[4>]), where E mn = E m 5 mn is diagonal and P[4>] a polynomial in (p, 
satisfies the system of Ward identities 

f d 2 Z dZ dZ \ 

We shall make a technical 

Assumption 2.2 The map i : I — > R defined by t(m) = E m is injective. 

This assumption is in fact less restrictive than it seems. For fixed potential tr(P[0]), 
the partition function is completely determined by the spectrum cr(E) of E and the 
spectral density. Since E is selfadjoint and has compact resolvent, cr(E) is discrete and 
has finite-dimensional eigenspaces. Moreover, by achieving E to be diagonal we have 
already placed ourselves in an eigenbasis of E. The requirements on E thus imply that / 
is necessarily countable, and the chosen diagonalisation yields a bijection between I and 
the set {(e, n £ ) : e G cr{E) , n £ e {1, 2, . . . , dim(ker(i? — e))}. This means that index 
summations over m G I of a function / which depends only on the spectrum of E but 
not on the multiplicity can be partitioned into 

dim ker(E— e) 

J2f( m )= E E f^ n e)= Yl dim(ker(£ - s)) ■ /(e) . 

mel eecr(E) n e =l e£cr(E) 

Now the new summation over o~(E) 3 e with measure dim(ker(i? — e)) is automatically 
injective. In other words, the Assumption |2]^ already covers the general case if we include 



an appropriate measure in the index summation. In order to make use of Proposition ^TT 
we need at least two indices a ^ p with E a ^ E v . This means that the only case which we 
really exclude is the standard framework of matrix models tr(E<ft 2 ) = /i 2 tr(</> 2 ) for which 



there are highly developed technologies |DGZ95|| . Note that this E = fi 2 idn does not have 



compact resolvent on infinite-dimensional Hilbert spaces. 
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2.3 Expansion into boundary components 

Feynman graphs in matrix models are ribbon graphs drawn on a genus-g Riemann surface 
with B boundary components (or punctures/holes/marked points/faces) to which the 
external legs are attached [LrWU5a|. Diagonality of the external matrix E mn = E m S mn 
guarantees that the matrix index is constant on every face (i.e. single line of the ribbon). 
This means that the right index b of a source J & coincides with the left index of another 
source J& c , or of the same source Jbb- Accordingly, there is a decomposition of W[J] into 
products of J-cycles J pip2 J P2P3 4 4 4 Jp n -ip n Jp n pi =: Jp n i where P n — Pi ■ ■ -p n stands for a 
collection of n indices. The decomposition according to the longest cycle reads 

oo oo N j N n j 

w[j] = w[o] + e e (n ^-) e ^i^m^i n n ^ ( 2 - 9 ) 

JV=1 ni,...,njv=0 j=l 3 J pi.eli j=lij=l 3 

The total number of J-cycles in a function G\ P i\_\ P i \ is its number B = 

ri\ + ■ ■ ■ + rijv of boundary components. For convenience we also list those terms of 
this decomposition which contain at most four sources: 

W[J] = W[0] + E G H t/ pp+ E {2 G \pi\ J pi J ip+ 2 G \p\i\ J pp J w\ 
pel pa^i 

~\~ ^ ^ ^ ~ G\p\q\r\ JppJ qqJrr ~^G\p\q T \JppJq r J r q ~\~ ~G\p qr \Jp q Jq r J r p ^ 
p,q,rEl 

^ ^ ^ ~^^G\p\q\ r \ s \JppJqqJ rr J ss -\- ~^J\p\q\rs\JppJqqJrsJsr- ~^G\ p qj rs \J p qJq p J rs J sr 



p,q,r,sel 

~^ G \ 1 P\'l rs \'^PP'^'l r '^ rs '^ S Q ~^*\p<irs\JpqJ qr-JrsJ sp ^ 0(J ) • (^'^^) 



In the sequel it is important to understand that the decomposition into cycles also 
persists for coinciding indices. For example, G\ pp \JppJpp is topologically different from 
G\p\p\Jp P Jpp. We keep track of this distinction by a continuity argument. For that we use 
an appropriate embedding I C M. k and a corresponding extension of E m to any C l -f unction 
(continuously different iable) on M. k with local extrema outside /. Perturbatively, the 
coefficients G\\ are rational functions of E m so that they become differentiate functions 
of the indices, too. Identifying the connected functions by functional derivative of W with 
respect to J a b it is then important to perform the derivative at generic indices, e.g. 



n d 2 W 

G\ pq \ 



(97 (97 

uupquu qp 



d 2 W 

J=o,q¥=p ' dJppdJ qq 



J=0,qjtp 



The case of coinciding indices is then obtained by continuity, G\ pp \ = lim^p G\ pq \ and 
G \p\p\ = lim g -> p G|p| g |. 

We are now going to prove the following Ward identity, which is the key to solve our 
model. 
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Theorem 2.3 For injective m i-> E m , the partition function Z[J\ of a real scalar matrix 
model with action S[4>] = tr(E<fi 2 + V[<j>]), where E mn = E m 5 mn is diagonal and V[4>] the 
trace of a polynomial in <fi, satisfies the system of Ward identities 

,_ r VJan'-JJnp • C/ p - C/ a _ r v UJ an UJ np / 



N N n j 



»?m-e e (ni) e nil- 

W=Oni njf=0 3=1 J J pi e p j=l £,-=1 



(G| a |a[pi|.„[P^| + ^2 G \Pl\...\Pi 1 \an\Pl\...\P^ N \ 
nel 

N 

+ Y1 Yl G \Pl\...\P^ 1 \nan qi ... qk . 3 \P^\...\P^ N \ J nqi-Jqig2 • • • J « 



fc=3 n,q±,...,q k _-j,el 
oo oo N 



<^>u j* lit, -• 

"M- e e (n-^)(n-J 



nAj n i ) V J - J - mdk mk 

N,M=0 ni,...,n N ,mi,...,Tn M =0 3=1 J fc=l 
N rij N m k 

yi Yi (nn ; / )(nn ^ ) ^ i - 1 ^ i - 1 ^ i ^ i - 1 « ^ ■ • • i « ^ i } • 

P? Qf <=I k J=l«i=l J fc=lZ fe =l 

Proo/. By definition, ]T -^LL = £ (^-^- + 1 > [4 We write 



r) f r) f p a L 1 ap 11 r) T F) J p a 1 J a P 

u J an u J rip tJ,J an u,J np 

where 5 ap VK*[J] contain all terms in which the derivatives generate a factor 5 ap , and 
Wa P re9 [J] are the remainders. Applying the functional derivatives to ( |2.9D we derive the 

structure of W^, W% as given in the Theorem. Namely, a 5 ap arises in Q j w ^ if (in that 
order) 

• both Try— and 75-7— act on different J-cycles each of length 1, or 

OJan OJnp 

• both Trf— and -Jj— act on the same J-cycle of length 2, or 

OJan OJnp 

• ttj— acting on a J-cycle of length k > 3 gives the chain Jpq J qoqi ■ ■ ■ Jq k _ 3 n- The 
second derivative -J^— then acts on J„ 0n in that chain. 

In ^y^ ^y^ , the only possibility to generate a 5 ap is to let 

OJan OJnp 

• both Try— and Try— act on J-cycles each of length 1. 

OJan OJnp J & 

The combinatorial factors in ( |2.9| ) are reproduced after a shift. It is straightforward to 
also write down the remainders W^ p re9 [J] in explicit form; but as we do not need this we 
refrain from listing the result. 
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We consider the case a ^ p. On one hand this implies S ap W^ = 0, on the other hand 
we may divide the expression in Proposition [D] by E p — E a ^ to obtain 



=> De + O^EAt^-il 3 )' (2-12) 

r ^ r - C/ p - C/ a v '-''Jan UJnp / 



It is now important to realise that the identity ( |2.12[ ) also extends to p — a. Namely, the 



partition function and its derivatives such as W^ re9 , W^ eg are completely determined by 
the entries E m in the external matrix. Extending m h- > E m to a differentiable function, 
the rhs of (|2.12|) is at least continuous at p — > a by l'Hopital's rule and agrees with the 



lhs for which there is no obstruction to go to p = a. This finishes the proof. □ 



3 Schwinger-Dyson equations 

3. 1 Strategy 

A standard procedure in quantum field theory consists in expressing the fields in the 
potential V[<p] = tr(P[0]) of the partition function ( |2.2| ) as derivatives with respect to the 
sources, 

d 1\ 

exp ( - ti(Ecp 2 ) + tr(0J)) . 



Z[J] 



V[<t>] exp - V 



Id J 



In this formulation the GauBian ^-integration can formally be carried out and gives 

Z[J] = C7e-^] c i«J>- , (J, J) E := £ . (3.1) 



The constant (possibly ill-defined) prefactor C will be omitted. The power-series expan- 



sion of e 



gives rise to Feynman graphs, but since the expansion looses the distinction 



between e v ^&j\ and e +v ^ej\ the expansion does not converge. A better strategy is to 

act on Z in order to 



'sjI intact and instead let functional derivatives 



keep e — u ^ """^ luul " UUJlu ' UW " U,,1 ' UJ vpg — aj 9p 

produce the connected functions G in the expansion ( p.9| ). These "external" derivatives 
combine with the "internal" derivatives in e~ v ^^ to certain identities called Schwinger- 
Dyson equations. It turns out that for each function G the Ward identity of Theorem [T3] 
can be used at an intermediate step to generate new relations. 

We shall demonstrate this for the regular (i.e. B — 1) two-point function Gi^i. As 
usual we impose a ^ b so that = \ , „. We have 

d 2 (log Z)[J] 1 



J=0' 



G|ab| 



9 J bad J a b 
1 f 9 



a 2 z[j] 



Z[0}\dJ, 



j=o 



z[o] <9j 6a dj, 

d 



ab 



1 &2[J] 1 dZ[J] 



j=o 



z[o] dj ba j=oZ[o] dj ab 



J=0 



ba 



+ 



1 



>ab- 



J=0 

9(-y) 



E a + E b 1 (E a + E b )Z[0] 
For the most general form of the potential V[(j)\ 



'ab 



d_ 

dJ 



Z[J] 



(3.2) 



K 

E 

k=3 pi,...,p k el 



E 



"T" ( / ) piP2 



J pk-lpk^pkpl 



we 
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have 



>ab- 



d(-V) 

d(fiab 



d_ 
dl 



K 



E E 

fc=3 n,p,p 1 ,...,p k ^ 3 £l 



k-2 



d 2 



dJiyni • • • d Jpi i _^b^ Jba dJan^Jnp 



(3.3) 



which means that in any part we encounter the two-fold derivative 



d 2 



9J an dJn 



known from 



Theorem [2.3| . This appearance of 
G. 



is a general feature of any connected function 



3.2 Schwinger-Dyson equations for a (ft 4 -model: B = 1 cycle 

We now specify to the case V[(f>] = ^tr(0 4 ). We first evaluate (|]| using; (13. 31) and 
Theorem |2l|, taking a ^ b into account and the fact that G\ q \ = 0: 

r 1 , (-A 4 ) td\{Wl[J] + Wt[J])Z[J\) 

(j\ab\ — — — — + 



E a + E b (E a + E„)Z[Q] I dJ ab dJ ba 

x ^ 1 d 2 f T dZ dz 

+ 1^ K _ K F) J ,F) J, [ Jpn r)T Jl 



^ E p — E a dJ pb dJf, a \ dJ an dJ np J ) j=o 

+ G\ a \ a \ ab \ + G| an | a b| + G\ aaa b\ + C|6a,6a| 

| G\qb\ - gjpfcj | G\g\ b \ - G\ b \ b \ \ 

"~ E p — E a E b — E a J 

To obtain the last line one has to use -^f^ = Z[0] [G\ an \J na + 5 an Yjq^i G\a\q\J qq + 0(J 2 )) 

and it = m(G\pn\J P n + K E qe i G \P\7 J gg + °( J2 )) ■ 

In a similar manner we obtain for pairwise different indices a, hi the following equation 
for the regular (i.e. single-cycle) (N > 4)-point function: 

r ("A 4 ) f d N {{Wl[J] + W 2 a [J})Z[J]) 



(_E a + i? 6 )ir[0] I dJ abl dJ blb2 . . . dJ bN _ 2bN1 dJ bN _ ia 

+ 1^ r _ r p t , ft i , . . . b i ; ftj, \ Jpn ~f)i Jna 



^ E p — E a dJ pbl dJ bl b 2 ■ • ■ dJb N _ 2 b N _ 1 dJb N _ ia ^ dJ np J J j=o 

Since all indices a, foj are pairwise different, the first line gives W„[0] + W„[0] times the 
application of all derivatives to Z [J] plus Z[0] times the application of all derivatives 
to (W^[J] + W„ [J]). In the first case there is actually only a contribution from the ap- 
plication of all derivatives to log2[J], because the indices are pairwise different. The 
second line is a sum of terms where one of the derivatives acts on J pn minus the action 
of Jb N _ ia on Jna, with all other derivatives acting on and ^r^, respectively. De- 
pending on the Kronecker-<5 arising in this way, the derivatives acting on Z [J] either form 
a single cycle or two cycles. For the derivative — ^ — — ^ special care 

to J J dJ P b 1 oJb 1 b 2 -dJ bN _ 2 b N _ 1 dJ bN _ lP 
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is needed in the cases p = bk- Either these are the necessarily arising cases in derivat- 
ives which form a single cycle, or the derivatives form two cycles 777 — a f . and 

777 |t — — ? with one common index bh. The single-cycle cases reduce to its 

9Jb k b k+1 - 9Jb N-2 b N-l Jb N-l b k 

action on log Z [J]. For the case of two cycles there are two possibilities: Either both 
act on a two-cycle contribution of logZ[J], or on one-cycle contributions of ^(log Z[J]) 2 . 
With these remarks we find 



G\ab!...b N -l\ - E +E b i ( G|a|a| + G \o-n\) G \ab 1 ...b N . 1 \ + G\ a \ a 



|o6i...6jv-i| 



N-1 



+ ^ ] G\ an \ a b 1 ...b N _ 1 \ + G\ aaa b 1 ...b N _ 1 \ + ^ ] G\b 1 ...b k ab k ...b N _ 1 a\ 
n£l k=l 

f 1 f 1 N— 1 ^/ y^i 

E (j |pfei...fe JV _i| <- J |afti...fe JV _ 1 | _ V"^ ^\b 1 ...b k \b k+1 ...b N ^ 1 b k \ { ^\b 1 ...b k \b k+1 ...b N _ 1 a\ 
R - R 2-^ 



E r , — E n Eh, — E„ 

pel p a k=i bk a 

G\hm 1 1 ...h*.r 1 b2i\ ^ 
&2i - E a 



JV-2 
2 



En ( - T |f>2I + l---&iV-lM (jr |&2I + l---&JV-ia| I /q r\ 

G |fei...b 2i | p p f- O 5 ) 



In Appendix || we derive the Schwinger- Dyson equations for functions with B = 2 
boundary components. 

3.3 The genus expansion 

Perturbatively, the connected functions G\ p i\^ p n \ are expanded into ribbon graphs, the 
matrix analogue of Feynman graphs. A ribbon graph represents a two-dimensional sim- 
plicial complex and as such encodes the genus g of a Riemann surface on which the ribbon 
graph can be drawn. We thus decompose the connected functions according to the genus, 
G\pI\...\pn n1 -^ G ip1i ^ n1 . 

The Schwinger-Dyson equations of the previous subsection describe how a given con- 
nected function G\ a \ P i\ ^\ P N | or G^ abpi pk 2 |pfc+i| \ p n \ is represented as a contraction of the 

special (planar) vertex ( gl ) or (0 aPfc _ 2 • • • )> which carries the external source 

J aa or Jab, to a smaller subgraph. This contraction may increase the genus. We have 
investigated a similar problem in [|GW05a and use these results to identify the topology 



of the various contributions to the Schwinger-Dyson equations. We can summarise the 
outcome of ||GW05a|| as follows: 



Proposition 3.1 i) Attaching one leg of a (planar) vertex to a connected subgraph 
does not change the genus. If this attachment is to a J -cycle of length 1, then the 
other legs of the distinguished vertex form a cycle of its own. If the attachment is to 
a J -cycle of length > 2, then the other legs insert in cyclic order into the connected 
chain. 

ii) Joining two legs of different cycles of a connected graph increases the genus by 1. 
At the same time, these two cycles are connected to a single one, unless both are of 
length 1 in which case both cycles disappear. 
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Hi) Joining two legs of the same cycle of a connected graph does not change the genus. 
The cycle disappears if before the contraction it was of length 2. If the original cycle 
was of length k > 3 and neighboured legs are joint, then the remaining legs form a 
single cycle of length k — 2. If the original cycle was of length k > 4 and legs are 
joint which are not neighbours, then the joining splits the original cycle into two, 
one consisting of the 'inner' legs relative to the contraction, the other one of the 
'outer 7 legs. □ 

We now identify the topology of individual parts of Ward identity and Schwinger- 
Dyson equations. We start with the two parts W„[J], W^fJ] in Theorem pT3| : W^[J] 
arises by insertion of {4> 2 ) pa into a connected function of W, and W„ [J] arises by joining a 
pair of disjoint connected functions from W ■ W via the vertex ((f) 2 ) pa . In the application 
to Schwinger-Dyson equation later on it will be important that {<fi 2 ) pa is part of a larger 
vertex (0 2 ) P a(0 <Xir )(f> 4 ) similar to ( |3.3| ). The r legs of (0® r )(f>i) have to be taken into account 
in the cycle structure. We find: 

1. G\ a \a\Pl\...\P N I ar i ses by contracting two cycles of the same graph, each of length 
1, via (4> 2 ) P a- We are in the situation of Proposition |3J].ii) so that the genus first 



increases by 1. But if (4> 2 ) pa is actually part of a larger vertex (0 2 ) pa (0 r ) a b, then 
the remainder (4> r ) a b forms a cycle of its own although notationally not visible in 
G| a | a |pi| ...\pN |- If these derivatives in the cycle ((fi r ) a b contract to h of the cycles 
iPj 1 1, ... , |P^J, then the genus further increases by h. The total increase of the 
genus is thus by h + 1 . 

2. Gpi| | P i \ an \P 2 \...\p^ | arises by joining both legs of a cycle of length 2 via (4> 2 ) pa - 
According to Proposition |3.1[ iii), the genus first is unchanged. The original J- 
cycle disappears in this way, but the remaining part (</>® r )(f>i) °f the vertex forms 
again a cycle of its own. If these derivatives in (<fi® r )(bi) contract to h of the cycles 
|P]|, . . . , |P^|, then the genus further increases by h. 

3. Gpii ipfc-2 ,„„„„ „ ipfci ipjv i also has unchanged genus according to Proposi- 
tion |3~T|.iii) after connection of the cycle of length k > 3 via (4> 2 ) pa - However, 
the final cycle structure is more involved: Assume the vertex is (j> pn (l>na( ( j>® r ) (&») an d 
to be inserted into the cycle J qk _ iqk Jq k qiJqiq 2 ' ' ' Jqh-2ih-i- We first contract (p pn with 
Jqk-2ik-ii giving a finq k _ 2 8pq k -i- According to Proposition |3~T|.i), the legs of the vertex 
insert into the cycle, giving 

In the second step (p na is contracted with J pqk , giving a 5 ap 5 nqk . Since k > 3 and 
r > 1, the two legs <f) na and J pqk are not neighbours, so that according to Propos- 



ition p3|.iii) the cycle (*) is split into two cycles J n q 1 Jq 1 q 2 " ' Jq k - 3 n and (<j> )(bi)- 
When Gpi, , p fc-a ,„„„„ „ ipfci i p jv i is used in a Schwinger-Dyson equation, the J- 

^1 I ■■■Kn fc _2 l' la ' l yi---9fc-3Kl I ■■■Kn JV I 

derivatives encoded by (4>® r )(bi) necessarily connect to h > 1 different J-cycles, and 
this increases the genus by h. 



4. G| a | P i| | P at |G| a |Qi| ...\qm i arises by connecting two length- 1 cycles of different sub- 
graphs via {<p 2 ) pa - According to Proposition |3.1| .i), the genus is unchanged. The 
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remaining part (0 <g)r )(6 ! ) of the total vertex (0 2 ) pa (0 (X "')(f> 4 ) forms a cycle of its 
own. In Schwinger-Dyson equations, this cycle connects to h > 1 other cycles 
hi G\a\Pl\...\p^ N \ G \a\Q\\...\QM M \, thus increasing the genus by h. 

5- Jpn & qJ J ^ — J na g fj J ^ ■ Here we have to expand Z as polynomial in connected functions 
G and apply the above discussion case by case. 

For the single-cycle two-point Schwinger-Dyson equation fl3.4|) we thus have: 
r {g) _ S 9Q A 4 ( v-> r (g') r (g") , r (s , ) r 0i") 

\ ab \~ E +Eu (E +Ek)\ ^ U \ab\^\a\a\ + Z^ Z^ ^ |a&| 



\baba\ 
nel 



Similarly, the equation ( |3.5|) for the single-cycle (N > 4)-point function has the expansion 



°'|o6i...6 J v_ 1 | 

(-A 4 ) 



\ Z^ Ljr M Lr |abi...6 JV _i| ^ Z^ U '|a|o| U '|a6i...6 JV _i| ^\a\a\ah. 



E 4- Eh I ' ' ' an ' l ab i--- 6 iv-il / -> |a|o| |abi...6jv-i| |a[o[o6i...6jv_i| 

a 1 g'+g"=gnei g '+ g "+i= g 



N-l 

|on|oci...cjv-i| |aaobi...6jv-i| ~ 1 ~ / / \b k ab k b k+1 ...c N ^ 1 ab l ...b k _ 1 \ 

nel k=l 



E ^lpbi-feiv-il Ur |a6i...bj V _i| _ V"^ u '|6i-6fc|6 fc+ i-6iv-i6 fc | ^ \b 1 ...b k \b k+1 ...b N _ 1 a\ 
E„-E a ^ E bl - E a 

pel p a k=\ bk a 

G (9") -G^' 
1=1 g'+ 9 "=9 " 2i 

In particular, we have proved 

Proposition 3.2 The (unrenormalised) planar regular connected N-point functions 
G\abi b N x\ ^ n a rea ^ 4> 4 - ma trix model with injective external matrix E satisfy the system 
of equations 

i \ \ r {0) - r (0) 

^(0) _ 1 M \- r (0) M \- ^\pb\ ^\ab\ ,„ ox 

G W-E a + E b E a + E b G ^^ G M + E a + E b ^ e p — E a ' ^ 

pel pel y 

r (o) r (o) _fM 

M°) - \ V" |pbi-fcjv-i| LT bbi|'- T |afei...fejv-i| 

Lr |a6i...6 JV _i| ~ A4 Z^ E — E 

pel p a 

^ G {0) - G {0) 

, \ y> £(0) ^,(0) |&2i+i-&JV-ifc2il |6 2 i + i-6jy-ia| ,„ „x 

4 Z^ [6i---6 2I j rp _ rp ■ v ' ) 

1=1 &2 ' a 
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r (o) 



The equation (|3.9| ) arises from the g = case of Q3.7| ) by elimination of the term 
via Q3.8|) . The equation ( p.8|) is a remarkable closed equation for the planar connected 
regular two-point function alone. This was already established in [|GW09|| by a different 
method, and versions of fl3.8|) for self-dual noncommutative 4 -theory at a = b = have 
already been derived in |DGMR0~7|j . We have generalised the ideas of |DGMR0~T| to a large 
class of matrix models and obtained a self-consistent equation for gSJ, for any matrix 
indices a, b. Usual Schwinger- Dyson techniques give an open system of equations relating 
the 2-point function to the 4-point function, the 4-point function to the 6-point function, 
and so on. In our setting, the Ward identity allowed us to break this system up into one 
closed non-linear equation for the 2-point function alone and then an inductive system of 
linear inhomogeneous equations for the (N > 2)-point functions, where the inhomogeneity 
and the linear factor are known by induction. 



We stress that the two-point function is by definition symmetric, G^ w 
this is not manifest in 



(0) 
\ab\ 



G \bl\i although 



4 Self-dual noncommutative ^-theory 

4-1 Definition of the matrix model 

In order to improve the problems of four- dimensional quantum field theory it was sugges- 
ted to include "gravity effects" through deforming space-time. A simple example for such 
a deformation is the Moyal space, which arises by deformation of the Frechet algebra of 



Schwarz class functions by the action of the translation group [ Rie93 



{f*g)(x)= [ %Mf(x+±ek)g(x + y )S k »> , f,geS(R d ). (4.1) 

Here, = — 0* € M^(R) is a skew-symmetric matrix. Although not required by the 
general theory, for our purpose we need to be of full rank, which implies that d is even. 
The algebra A e = (<S(M d ),*) is a pre-C*-algebra ||GG1SV03| , Prop. 2.14]. The standard 



]R d -Lebesgue integral defines a linear functional which is positive L d (f * f)( x ) > an d 
tracial L d (f * g)(x) = f Rd (g * f)(x). Therefore, GNS-construction gives rise to a Hilbert 
space H = Hi g) H2 on which there are commuting representations it <g> id of Ae and 
id eg) n op of the opposite algebra Aq. Restricted to Hi, any element of Ae is a compact 
operator on Hi, in fact even a trace-class operator. 

This means that we can define field-theoretical matrix models on Ae 9 4> by specifying 
the polynomial P[4>] and the external matrix E. We take d = 4 and P[4>) = ^0 4 - The 
unbounded operator E should play the role of the Laplacian. However, the Laplacian A 
on R d does not have compact resolvent. This non-compactness of A is the reason why 
field theories based on action functionals such as 



/ dx (~<f){- A + n 2 )(p+ —0*0*0*0) 
J Rd \2 ' 4 / 



x 



have bad properties (UV/IR-mixing ||MVSOO|0 . 



In our previous work [ |GW05b|| we found a way to handle this problem. We realised 



that extending the Laplacian to the harmonic oscillator Hamiltonian (which has compact 
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resolvent!), the field theory defined by the action 



S 



-A + n 2 \\2e- l x\\ 2 + f?) 



X 

+ 4 ( 



(4.2) 



is renormalisable to all orders of perturbation theory. The model is covariant under the 
Langmann-Szabo duality transformation |LS02|| and becomes self-dual at Q = 1. Certain 
variants have also been treated, 



sec 



Riv07a for a review. 



Writing the action functional ( f4.2| ) in the matrix base of the Moyal space, see ||GW05b| . 
GW05a|| , one obtains a non-local matrix model in the sense that the 0-bilinear term is 
a tensor product ti(E ■ (</> <g> </>)). We have studied the power-counting behaviour of such 
matrix models in ||GW05a|| . By a sophisticated analysis of the tensor-product external 



matrix E we proved in [ GW(J5b| ] the perturbative renormalisability of the model 
The action simplifies enormously at the self-duality point Q = 1 where it becomes an 
ordinary matrix model in the sense of section |2| with diagonal external matrix E. In 
terms of the bare quantities, which involves the bare mass pbare and the wave function 

(2.5)- 



renormalisation d> i— >■ Z^ch, the action functional becomes GW05b 



eqs. 



-(2.6)] 



S[4>]= ^ E m0mn4 > nm + V 

E m = (2n9) 2 z(-\m 



t^bare 



V[</>] = (2n9) 



,Z 2 X 



(4.3a) 

^mn<firti<fikl(f>lm ■ (4.3b) 



We have made use of the fact that, after an appropriate coordinate transformation in M 4 , 
the only non-vanishing components of 9 are O i2 = — O21 = ©34 = — @43 =: Q- We have 
also absorbed the ground energy of the harmonic oscillator into a redefinition of fibare- All 
summation indices m,n, . . . belong to N 2 , with \m\ := m x + m 2 . The symbol Njj- refers 
to a cut-off in the matrix size which later on will be made more precise. The scalar field 

is real, 4>mn — 4>nm- 

Note that we use here the bare coupling constant A. In general we should be prepared 
to put A = Z4X ren for a suitable divergent Z4 and a finite X ren . However, according 
to | pGMK(]7[ the model under consideration has (perturbatively at any order in A) no 
infinite renormalisation of the coupling constant. In this paper we give a non-perturbative 
proof of that result. 



4-2 Renormalisation of the two-point function 

We now focus on the Schwinger- Dyson equation ( |3.8| ) for the planar regular connected 
two-point function G^,. From (|4.3b|) we read off A 4 = (2ir6) 2 Z 2 X. The summation 
over p G J := in ( |3.8|) diverges if we remove the cut-off N 2 ^ — > N 2 and, therefore, 
requires renormalisation. Renormalisation is defined by normalisation conditions for the 
first Taylor coefficients of the one-particle irreducible (1PI) function Tab related to 
according to 

G\% = 7^{H ab _ - r^)- 1 , H a _ b _ := = z(- e (\a\ + \b\) + p 2 bare ) . (4.4) 
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It is then an easy exercise to rewrite fl3.8| ) into the following equation for T^: 

XZ 2 ^ / 1 1 \ XZ ^ 1 r p6- r a6 



ab 



^ / 1 1 \ XZ ^ 1 

(2^)2 2^ U ap -r ap + f^-rj + (2^)* ^ (I^r p6 ) |(|p|-|o| 

- - - - E fc w 



(4.5) 

We require 

r,_ b = z»l re - /i 2 + (z - i)l (|a| + 16|) + rr , 

= (<9r en )oo = , (4.6) 
where 0T ren is any of the (discrete) derivatives with respect to ai, 02, &i, &2- This implies 

= tf a6 - T ab = ~(|a| + |6|) + /i 2 - . (4.7) 

Hence, \x is the renormalised mass, and both G[% and T r a e b n should be regular if the cut-off 
in the matrix indices is removed. Inserted into (kL5|) we find 



z^l re -ii 2 + (z-i)l(\a\ + \b\) + r 

A 



ren 
ab 



(2^) 2 ^|(|a| + b|)+^-r- 

A / z z ^ r P t n -Kb 



(2^) 2 ^m+\p\)+^-r;r ~ m+\p\)+^ 2 - T ?b n tl\p\-\a\)) ■ ( ' 8) 

Notice the different exponents of Z in the two tadpoles J2 P i^+^+^^ren ! Separating 
the first Taylor term (by putting \a\ = \b\ = 0) we obtain 

A ^ Z 2 



ran. 



A ( Z Z 



E 



(2t») 2 ^ V 1 |p| + ^ _ r ~» l\p\ + ^- rg» 7 • 



which inserted into ( |4.8|) gives 
(Z-l)}(|a| + |6|) + r3» 



A v / Z 2 Z 2 

(2^) 2 \ l(\a\ + M) + - r™ fbi + ^ - rg* 

A ^ / Z Z 



E 



„ pren pren „ pren 



+ IpI) + n 2 — r;r |(bl - Is!) |bl + a< 2 - I\p\ 
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The next step consists in differentiating (f4.9|) at with respect to a, and hi in order 



to get a self-consistent system of equations for Z, T r a l n . It is now crucial that ( }4.9|) 
depends only on the sums \q\ = ai + a 2 , \b\ = 61 + 62 and \p\ = pi + P2 of indices. 
This degeneracy is first passed on to the bare two-point function T ab . But then the 
derivatives £|y,„M = and & | |aU5M = 4l lsl ,, M respect the degeneracy, 

so that also the renormalised two-point function T r a l d depends only on the 1-norms 
\q\ = ai + a 2 andJ6| = b\ + 62- Consequently, we can replace the index sums by 
X^eN 2 f(\p\) — J2\p\=o(\p\ + l)/(bl)> because there are |p| + 1 possibilities^ to write 
|p| as a sum \p\ = pi + p 2 with pi,p 2 G N. The cut-off [p| < A/" G N then specifies the 
previously unclear symbol p G Njj-. In summary, ( |4.9|) takes the form 

(z-i)j(|a| + |6|) + rr 



A 



ab 



^ 2 1^0 v l^l + bl) + m 2 - r™ |b| + ^ - rgr- 

A ^ / Z Z 



E(H + i ; 



run 



„ pren pren „ p 

Z L pb L ab Z L 



ren 
pO 



+ iu . ... rrm in • ( 4 - 10 ) 



+ |p|) + /i 2 - r-« }(jp| - |a|) + ^ _ r -« §h 



5 Integral representation 

5.1 Wavefunction renormalisation 

We study a particular limit in which the self-dual noncommutative </>4-model converges 
to a large- M limit of a certain matrix model. This limit is defined by sending 6, M — > 00 
such that J 

-A/" = yu 2 (l + X)A 2 = const . (5.1) 
o 

Here 3^ > —1 is a real number which we identify later in order to simplify our equation. 
The number (1 + y) can be seen as a finite wavefunction renormalisation. This large- M 
limit (|5.1|) plays the role of a thermodynamic (infinite volume^]) limit because it turns the 
discrete norms \q\, \b\, \p\ of the matrix indices into real numbers a,b,p G [0,A 2 ] defined 
by 

}\a\=:f{l + y)a, ||6| =: /x 2 (l + y)b , f|p| =: ^(1 + y)p , =: /i 2 r afc . (5.2) 

In the very end we also have remove the cut-off A 2 in a continuum limit A — > 00. It 
is important that only depends on \q\ and |o| so that it converges to a function 
[0,A 2 ] 2 E5 (a, 6) r-> r afc G R. In the limit ( |5.1|) , a sum over |p| converges to 

A 2 



E /(f bi) *3 /i 2 (i + y) / * /(/(i + w • 

„l_n JO 



bl=o 

1-1 



1 In d dimensions the measure would be ( - d 2 1 



2c 



See [GW11, GW12] where it is shown that the spectral action behind the model under consideration 



yields for il = 1 an 8-dimensional finite volume proportional to 9 
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Therefore, the equation ( |4.10p converges to the following integral equation where the mass 
/j? drops out (or rather becomes 1): 



(z-i)(i + y)(a + b) + r ab 
x(i + y) 2 rA2 



pdp 



64tt 2 J \( a + p)(i + y) + i-r ap P (i + y) + 1 - r 0p 

A(l + ^) 2 f A ' ( Z Z 



64vr 2 J pdp {(b + p)(l + y) + l-T pb p (l + y) + l-T p0 



Z T pb — T ab Z T p Q 



(b + P )(i + y) + i-T pb (i + y)(p- a ) p(i + y) + i-r p0 P (i + y). 

(5.3) 

The powerful analytical tools available for integral equations are a huge advantage over the 
discrete equations. On the other hand, in regions of the space of continuous parameters 
where T varies slowly we can expect the integral equation to be a good approximation for 
the discrete equation ( 4.1C ) also for finite 9. 

It turns out that the terms with Z 2 -coefficient in ( |5.3|) need a subtle discussion. It is 
highly convenient to eliminate them via the equation resulting from ( |5.3|) at b = 0: 



(z-i){i + y) a + r a0 
\(i + y) 2 M " 



palp 



64tt 2 J \( a + p)(l + y) + l-T ap p(l + y) + l-T, 



Op' 



+ M±±yi[ A dp z aT p0 -pT a0 

The difference between (|5.3|) and (|5.4|), expressed in terms of 



Gab: ~ (a + b)(i + y) + i-r ab > (5 ' 5) 



^(1 + ^-^ + 7^ = ^ pdpZ^-^L. (5.6) 



reads 



G ab G a0 64vr 2 J p-a 

From ( |5.6| ) we obtain the desired equation for Z^ 1 by putting a = 0, dividing by Z(l + y)b 
and going to the limit b — > 0, where lim^o Kg^ — ^) = 1 + ^ an< ^ um 6->o G ;, = 1 are 
used: 

Z = 1 ra 2 / G p0 - im / dp . 

647T 2 J QAti 2 b^o J Q b 

It will be convenient to choose y as 

-A 2 r< r* \ /"A 2 



(5.7) 
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Of course it remains to verify that — 1 < y < oo. Multiplication of ( |5.(j| ) by b ^°y^ z and 
elimination of Z~ l via (|5.7|) gives 

(G- _ G-) (l + _L. ) + Ga0 = > r" dp ( p^-^-^-^gS ) 

V &GW 64n 2 b J \ p-a J 

(5.8) 

5.2 The Carleman equation 

We make now the crucial assumption that G ab is Holder-continuous, i.e. 

\G pb - G ab \ < C b \p - VO < a ± p < A 2 , (5.9) 

for some constants C b , r\ b > 0. Under this assumption we may replace the singular integrals 
by their Cauchy principal values: 



dp Gpb - Gab = lim f / +/ )dp g*^g* . (5.10) 



Introducing the finite Hilbert transform 



^[/(.)]:=-hm / +/ \^L, (5.11) 

K ^° V Jo Ja+e ' P-d 



we can rewrite the equation ( |5.8|) as 



Va aG a0 ' 647T 



Dab := a Gab r ° a0 • (5.12b) 



The equation ( |5.12a| ) is now a standard singular integral equation of Carleman type | JUar22 



Tri57|| . We cite from Tricomi's book |[Tri57|| : 



Proposition 5.1 ( || , l'ri57| , §4.4]) Let h E CQ-1, 1[) and f E L p (}-1, 1[) for some p > 1. 
T/ien i/ie singular integral equation 



h{x)<t>{x) - XirHaM*)] = f(x) , x e ]-l, 1[ , (5.13) 

/ias i/ie solution 

0( x ) = sin (^)) f cos( ^( x) ) + e «.[^ t [ e -«.W/(.) sin(^(»))] + ^1!!) , ( 5 .14) 

A7T ^ 1 — X / 

n/ n +■ f Xir \ . ... I Avrl ... sin($(x)) 

[o,-] ^(h(x)) 2 + (Att) 2 tan^(x)) 

where the Hilbert transform integrates over ] — 1, 1[ ; and C is an arbitrary constant. 
The angle d(x) obeys the identities ||Tri57| , §4.4(28)] and ||Tri57| , §4.4(18)], 



e -H x [#) cos(t?(x)) + n x [e~ n 'W s in(tf (•)] = 1 , (5.15a) 
e H.[*\ cos (>d(x)) - n x [e n ' w sin(tf(«))] = 1 . (5.15b) 
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The solution ( |5.14p and the identities ( [5.1 5a ) and ( |5.15b| ) are unchanged if transformed 
via p = 4r(l + x) to p G [0, A 2 ] instead of x G [—1, 1]. We make the following decisive 

Assumption 5.2 C = . 

This assumption (or rather choice) will be discussed at the end of Sec. 5^3 and in Sec. [7]. 
Under this assumption the solution of ( j5.12a| ) is 

D ab = -^sin(^(a))(G a0 cos(^(a)) + e^^^ [e~^ M G.o sin(^(.))]) , (5.16a) 



Mo) 



arctan 

[0,7T] 



1 + bG a0 + ^Hi [G.o] 



(5.16b) 



The form ( |5.16a| ) of the solution of ( p.!2a ) is not very useful. We can simplify it enorm- 
ously noting that ( |5.16b|) is, for 6 = 0, also a Carleman-type singular integral equation 



A cottf (a)G a0 --^^[G 



647T 647T 

This equation has the following solution: 



o0 



(5.17) 



^o1-^o1 sin (^o( Q )) 

|A|a 



Lemma 5.3 G a0 = QAne 

Proof. By ( |5l4j) the solution of ( gT7|) is for C = 
647r . , „ , ss /cos(-#o( a )) 



G 



«0 



A 



sin($ (a)) 



-MiPo] 



+ e « A a [^A ( "' ^sin^oH) 



Rational fraction expansion H a 



[/(•)] -^o [/(•)]) and the identity 
( |5.15a| ) yield the assertion, where ^o(O) = for A > and $o(0) = 7r for A < 0, hence 
cos^o(O) = sign(A), are used. □ 

For the Hilbert transform arising in (|5.16a|) and for the investigation of the four-point 
function later on we need the following addition theorem: 

Lemma 5.4 For all < a, b, d < A 2 one has 

Xa 



64tt 



sin 



v&d{a) — $b{a)) — (b — d) sin$&(a) sint^a) 



(5.1* 



□ 



Proof. This follows from insertion of ( p.!6b| ) into coti9{,(a) — cot^(a). 
This Lemma has several important consequences: 

Corollary 5.5 1. For X > 0, the function b i— > $&(a) is monotonously decreasing for 
any fixed a > 0. In particular, $&(a) < 'do (a) f or an U < a, 6 < A 2 . 

2. For A < 0, the function b (-> i) b (a) is monotonously increasing for any fixed a > 0. 
In particular, i?&(a) > $o( a ) f or an U < a, 6 < A 2 . 
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Proof. By continuity we can choose \b — d\ small enough so that r a M : = ^(a) — $&(a) G 
] — |, | [. Then r a ^ and sinr &d have the same sign. Let b > d. Then from (|5.18|) we 
conclude i? d > $ b for A > and < i? 6 for A < 0. □ 

The next step is the analogue of Lemma pT3| for the Hilbert transform arising in ( |5.16a ): 
Lemma 5.6 

H A [ e -u\[M G%Q S in(^(.))] = ^®e"* ( e ^[*°-*] cos (tf ( a )_tf 6 ( a )) _ 1 
Proof. We insert G.o from Lemma |5.3| into 7^ [e~^*^G.o sin($6(«))] and obtain with 



sign(A) 



sin (tf (.)-tf 6 (.)) 



Now the the identity ( |5.15b[ ) yields the assertion. 

This allows us to prove the following remarkable formula for G a t,: 



□ 



Theorem 5.7 The renormalised planar connected two-point function G a b of self-dual 
noncommutative (pf-theory (with continuous indices arising as limit 9 — > oo) satisfies 
( under the assumption \5.3j) the equation 



G„h = 64vr 



sin(^(a)) M h] ^ m 



A a 



(5.19) 



In particular, G a b > for all a, b G [0, A 2 ]. 
Proof. Insertion of ( |5.12b| ) into ( |5.16a| ) gives 

MM*)) ( Ga0 cos(Ma)) + e n*M H A e -^[^I G . sHM9)) 



Gab = G a0 - 647T6- 



Aa 



We insert Lemma 5.3 and Lemma 6.61: 



Gnh = 647re 



■^W^ife*] - 64nb Sin{Ma)) \e^ cos(^(a)) 



A a 



Aa 



+ e 



Aa 



64nb sin($o(a)) 



Ml#o-# b ] cos (0 o ( a )-0 6 (a)) - l) 



We express cos (^o( a ) — $b( a )) by its addition theorem and combine the sin-sin part with 
e K\M a f ter t he first "{": 



G nh = 647re 



-«A W ?MM( e «iN cos 2 (^(a)) + c ^f^l sin (^( Q )) 



A a 



sin(# (a)) 



e Km sin (tf 6 ( a )) cos(^ b (a)) f ^ + cot o (a) 



Aa 



Now ( B.16b|) implies _|_ coti9 (a) = cot-# b (a), and then the above equation collapses 
to " □ 
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Lemma 5.8 One has 




(5.20) 



Proof. To compute y 



A 



A 2 



dp 

647T 2 J p 



D 



A 



pf) 



64vr 



T-Lq[D, ] in (|5.7|) we use the convolu- 



tion theorem for the finite Hilbert transform |[Tri57| , §4.3(4)] 



»)n.[<h] + M*W 



With D p o given by (|5.16a| ), we have to identify x — 0, </>i(a 
02(a) = e-^°]G a0 sin(??o(a)). With sin$ (0) = we obtain 



?4[<M»)]?4[0 2 («)] - 2 (x)0i(x) . (5.21) 

: e ^a[^o] S i n ^ ( a ) a nd 



n 



+ ^[ S in^ (.))e w -^]^[ e -^^G.osin^o(«)) 
Using (|5.15b|) this equation collapses with cosi9;,(0) = sign(A) to 

y = sign(A) \e- HAM+Ho[ » o] G. sin(tf («)) 

• A '-' sin 2 (^ (p)) 



sin(A))e w -H e-^^G.o sin(tf («)) 



64 
T 



dp 



dp 



6 ^ 2 J° (^) 2 (l + cot 2 (^ (p)) 



where the second line follows from Lemma |5.3j . Now the assertion follows with (|5.16b|) 



□ 

It follows that y > for A > and — l<3^<0forA<0 with |A| small enough. Now 
observe that fl5.7| ) , Lemma |0] and ( |5.15b| ) give 



^-^(l + ^^-sig^A)- 



^olM r A2 dr> 

— / sin^ (p)e^^ : 

K Jo P 



l+y)(l- sign(A)e-^^o] H A [sin?9o(#)e ^M]\ = ( 1 + 37) si g n (A) e 



-n A [M _ 

(5.22) 

It therefore follows that sign(Z _1 ) = sign(A) so that the case A < leads to an unstable 
action functional. We therefore discard the case A < in the sequel. 

5.3 The master equation 



The identity (|5.19|) allows us to compute G a b once G a0 is known, but ( |5.19|) alone is not 
sufficient to determine G a Q. This is because we ignored so far the equation (|5.4|) which 
in terms of (15.51) and after insertion of (15.71) reads 



G 



aJO 



\{i + y)z 

64tt 2 



pdp (G ap - G 



A 



op; 



64tt 2 



dp 



a — a 



G 



GaO 



(p-a) 



(5.23) 
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We introduce the Hilbert transform and insert ( 5.19p and (|5.22[ ): 



1 + a - 



Xa 
64^ 



cot $0(0) = — pdp 



A 2 



sin $Ja) 



sin fiJa) 



7ra 



Aa 
64vr 



(5.24) 



We write sin$ p (a) = (1 + cot 2 $ p {a)Y^ = (l + + cot $o(a)) 2 ) 2 and use the following 
limit for A > 0: 



lim r ®n(M a )) 

a~^o V na 



cos-i^O) lim 



P (a) 



A 



a^O a 647T 2 (1 + p) 

This gives: 

Proposition 5.9 The function T a '■= ^r c °t^o( a ); w ^h % — 1, is determined by 



(5.25) 



T a = 1 + a 



A 



n I! 111 

6471-2 e _> 



a+e 



dp 



a 



P exp H A a 



p — a 



+ 



arctan 



P+T. 



p exp H 



arctan 



P+T. 



l+p 



. (5.26) 
□ 



This is a non-linear self-consistency equation whose solution yields the input data $o(a) 
and $b (a) to evaluate the two-point function G a b by Theorem |5T7| . 

Unfortunately, ( |5.26|) is too complicated for further analysis. We therefore derive an 
alternative equation from the symmetry G w = lim a _> G a b- This requirement together 
with Theorem |5]7| gives: 

Proposition 5.10 The renormalised planar regular two-point function Gbo is determined 
by the self- consistency equation 



1 



A 



60 



1 + 6 



exp - 



64tt 2 J 



A 2 



dt I dp 



(GpoY 



which is well-defined in the class of continuously differentiable functions on [0,A 2 ]. 
Proof. The limit a in ( pTTSD yields with (|5~25D 



, (5.27) 



Gbo = lim Gab 

a— >0 



1 + 6 



1 + 6 



exp 



dp 
^ Jo P 



{Mp)-Mp))) • (5-28) 



Together with Corollary |5.5| this formula already establishes that for A > the function 
b (->■ Gbo is monotonously decreasing and maps [0, A 2 ] into [0,1]. The explicit formula 
fl5.27T) follows from 



- f — («o(p) - Aip)) = - f" — f dt 

^ Jo P ^ Jo P Jo 



dMp) 
dt 
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and (5.16b) after exchange of the p, t-integrals. 



The equation ( |5.27| ) is meaningful on the class of continuously differentiable functions. 
Namely, if G m0 G C^O, A 2 ], then 



64vr 

\p 
64^2 



^a [g . o] = _^ , Lmo ( f + y q Gj0 | io g i, - p| 

(^G A 2 log(A 2 - p) - log p - dq G' q0 log \q 



P\ 



This shows that p h-» ^"HpfG.o] is a continuous function on [0,A 2 [. Since 
lim p _ s>A 2 ^'HpfG.o] = — oo and ■^'H^ ) [G. ] appears in the denominator together with 
bounded functions, the integrand in ( |5.27| ) extends to a continuous function on [0,6] x 
[0, A 2 ] 3 (t,p). By the fundamental theorem of calculus we have 

- -G» (tA + A, ^ , . „ - , i ) ■ (5-29) 



a* VI + 6 64** Jo (^ Gjj0 ) 2 + (i+6Gjj0 + ^a [g . o]) - 

Since the rhs is a smooth function of 6 G [0,A 2 ], continuous derivatives of G(, exists 
inductively to any order. □ 



Observe from (|5~29|) and ( pOj ) that ^| fe=0 = -(1 + 3>). 



The next task is to prove that the master equation (|5.27|) has a solution for any A > 
and A 2 sufficently large, including the continuum limit A — > oo. This will be achieved by 
the Schauder fixed point theorem: 

Proposition 5.11 (Schauder) Let X be a Banach space, K C X a convex subset and 
T a continuous mapping of K into itself. IfT(K) is contained in a compact subset of K, 
then T has a fixed point. 

The main road to establish compactness for X being the space of continuous functions 
on a compact Hausdorff space is the Arzela-Ascoli theorem. There exist generalisations 
to differentiable functions and to locally compact Hausdorff spaces. The setup relevant 
for our situation is contained in |(J(JMX1"T[ | that we specify to our case. Let 



C (M + ) := {/ : M + — >■ K continuous and bounded , lim f{t) = 0} 

t— >oo 

C 1 (M+) := {/ : K+ R : /, /' continuous and bounded , lim f(t) = lim f'(t) =0} 

be the spaces of continuous and continuously differentiable functions on R + , respectively, 
that vanish (together with their derivative) at oo. These are Banach spaces with norms 
H/lloo = sup 4eR+ \ f{t)\ and WfW'^ = sup tm+ \f(t)\ + sup teM+ |/ ; (t)|, respectively. One has 

Proposition 5.12 ( |[Avr69| , Lemme 1]) A subset K C Co(M+) is relatively compact if 
and only if it is uniformly bounded, equicontinuous on R + and uniformly convergent at 
oo. 

We recall that a subset K C X is relatively compact in the topology induced from X if 
its closure K is compact. A subset K of continuous functions on a metric space (M, d) is 
equicontinuous if for every e > there is a uniform 5 > such that \ f(t) — f(s)\ < e for 
all / G K and all s,t G M with d(s, t) < 5. For subsets of Cq(R + ) one has: 
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Proposition 5.13 ( JGUMXTX] , Thm. 3.1]) A subset K C C 1 (M+) is relatively com- 
pact if and only if 

1. K' :— {/' : / G if} is relatively compact in Co(M+), 

2. For every e > t/iere is a L > wift ||/| j L J]^ < e /or a// f E K. 

The following operator T A is well-defined on the subset of positive functions in Cq (R + )D 
L q (R+), for any 1 < q < oo: 



m/)(6) : . _1_ exp ( _ _L ^ f _ . - » - ] . (5.30) 



The Hilbert transform maps differentiable functions to continuous functions on ]0, oo[, 
with logarithmic divergence at p — 0. To see this we write for p > 

WW lJ/%MiiM + ir,,M, (5 . 31 ) 



using the fact that T-L^[l] = 0. The second integral is bounded by Holder's inequality, 
and its derivative with respect to p exists. In the first integral we have 

2(P+<I dq /(,) " /(P + S) f " /(,) " /(P) 



*"> dq m- fi + j) I p dq I \ m-m + /(p + *)-/ M) . 

2p g-(p + d) y v g-(p + ^) g-p g-(p + o) / 

Using the mean value theorem f(p) — f(q) = (p — q)f'(£) for some £ between p and g, this 
expression tends to zero for 5 — » 0. Since linip^oP^p 3 [/(•)] — 0, the integrand in ( |5.30| ) 
is a rational function of positive continuous functions on M + with denominator separated 
from zero. Therefore the p-integral exists and is convergent for large p. The fundamental 
theorem of calculus guarantees that (T\f) is differentiable. 
The following even stronger result holds: 

Lemma 5.14 Let A > and assume that f G Cq(M+) satisfies 

(Al) /(0) = 1, 

(A2) < fib) < jL., 

(A3) 0<-(^ + £f)<(7 A , 

for some C\ > 0. Let P\ > be the unique solution of + C\)e° xPx = 1. TTien 

(CI) (T A /)(0) = 1, 
(C2) 0<(T A /)(6)< T ^ ; 

l^Oj U S U+l + (T A /)(6) J - 2 + \P X > 
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(CA\ | ( r A/)"( fe ) | < 23 | 2 | 7+87T / 64 \2 
l^J I (T A /)(6) I - 4 "T tt "T 2 \\PJ ■ 

Proof. (CI) and (C2) are obvious from ( |5.30| ). 

Integrating (A3) with initial condition (Al) one has e 1+ ^ < /(6) < which is even 
stronger than (A2). In analogy to ( |5.29|) we have 



(T/)(6) 1 + 6/ 647r 2 7 ( 6 A / ) 2 | (6 | 1+ ^^')] ) 2 ' 

The rhs of (|5.32| ) is positive, and it remains to show the bound (C3). Since / > and 
/' < we have the following inequalities for the integrals in ( |5.31| ): 

< / dq < — — log , 

hp Q-P P+ 1 P 

> f P d q f ^ - > 2p inf /'(g) > -2p(C x + 1) . 

JO q-p 96[0,2p] 

Collecting these results and using s jt^~ < /(&) < jt^ we conclude 

1 + > (1 - ^(1 + ^)e^) . (5.33a) 

For fixed Ca > the continuous function M + 9 p i— > ^?(1 + C\)e CxP is monotonously 

increasing and maps M + to itself. We let P A be the unique solution of ^^(l+C\)e CxPx = 1. 
By the implicit function theorem, this gives rise to a globally defined different iable function 
R + 3 C\ h- > Pa which is monotonously decreasing. These monotonicity properties lead to 
the inequality 

<P< P X =>. L±fcW(!)] >i_^>q, ( 5.33 b ) 



/(P) ~ ^ 



We now split the integral (|5.32 ) at P\ where for < p < P\ we use ( |5.33b| ) whereas 
for p > P x we simply use (b + 1+ ^| [/(,)1 ) 2 > 0: 



(T/)(6) l + 6/-64^y (^) 2 +(l-^) 2 64^ y Q 

AP A Z" 1 ete 64 



64tt2 7 ^)2 X 2 + {1 _ X 2 ) 2 + X P x 

\Px_ f 1 dx(l + x 2 ) _6£ = 1 _64_ 

" 647r2 /o (^) 2 x 2 + (l-x 2 ) 2 AP A 2 + AP A ' l ' J 



Z 1 1 + x 2 1 Z" 00 cosh I 

We have used the integral / dx — r— — ^— = - dt 



o 



TT 



a 2 x 2 + (1 - x 2 ) 2 2 A f^-n + cosht 2a 



which can be looked up in |GR94j , §3. 517.1, §8. 702]. This finishes the proof of (C3). 
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Another derivative of fl5.32| ) gives 

1 X f°° dp 



(T x f)"(b) = (Tf)(b) 



1 + 6 ' 64vr 2 7 (Ap) 2 | (ft | ^±jfcWy 



1 « 2 (ft + l±jfcg[/(!)h 
+ 7X1 + «TT / * 7 — /( f j ^ - (5-35) 

We thus get with (C3) the estimate 

|(T A /)"(6)| < /3 + _6£\ 2 + 1+ 2A r dp 



(T x f)(b) ~\2 AP A V 647r 2 i ^ Ap , 2 | ^ | 

/3 64 \ 2 / 64 \ 2 XP X f 1 dx 

< L + ttH +1 + 4tt — - + 



2 AP A / VAPa/ 32tt 2 



o ((M^ 2 + a-* 2 ) 2 ) 



We have with ||GR94i §3.517.1, §8.702, §9.131.1, §9.122.1] 

f 1 (l+x 4 )dx _ 1 r°° cosht 
Jo ( a 2 x 2 + (l-x 2 ) 2 )~ 2 v^io ((f-l)+cosh*)l 

"4V 2Fl l 2 



3tt ^\3jt pf-is 

4 / 8a 2 1 V 2 



l--,K 1 + a 



a 2 / a 2 



Together with 2a < 1 + a 2 we confirm (C4). □ 

It is < 
Namely, 

^ + AP7" A W 32^V2"apT; C " - A "32^' 



It is of fundamental importance that we can always choose C\ such that | + < C A . 
1 64 „ AP A 2 /3 64 \ ( i+j^)p a AP A 2 ,„ „ , r p 



* Pa6X - w + ^) e " <1 - (5 - 36a) 

The condition is satisfied for 



(14 



P, = ^= =» P A e¥.(^ + |)e*< 0.35714 <l. (5.36b) 



Note that the inverse solution for C\ is huge for < A «C 1. 
We thus consider the following subset of Cq(M + ): 



{feC 1 (R + ):f(0) = l, 0</(6)<-^, < -f{b) < + C> 



(5.37) 



AP? 

)iution or 

following facts 



. I>4 



with C\ the solution of ^^-(l + C\)e CxPx — 1 at P\ — -, x A . We have established the 

V 1+ 64 
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(Kl) T\ maps K\ into K\ by Lemma §g(Cl) + (C2) + (C3) and 

(K2) (T A (/C A )) C C (M+) is relatively compact. Namely, (T A (/C A )) C C (M+) is equicon- 
tinuous by Lemma |5.14] .(C4), uniformly bounded as subset of /C A and uniformly 



convergent by Lemma |5.14[ (C3) + (C2), thus relatively compact by Proposition |5.12| . 



(K3) \\gL \\ < e for all g G T X {K X ) by Lemma CT(C2). 



(K4) For the closure we have T A (/C A ) C /C A : Functions g = lim^oo <7fc with g>& G T A (/C A ) 
on the closure can at most exceed /C A in g(b) = for some b. But < —g'{b) < 
(jj^ + C\)g{b) and g(0) = 1 imply e 1+ ^ < g(b) < •jtj so that also g is strictly 
positive. 

(K5) fC\ is convex: We have /i/i + (1 — /x)/ 2 G /C A for any / 1; / 2 G /C A and < /i < 1 by 
definition of /C A . 

(K6) T\ : /C A — >■ /C A is continuous, i.e. ||T A / - TJW^ < e for ||/ - fW^ < 5: By the 
uniform convergence Lemma |5TT4l. ( C2 ) + ( C3) with \ + jp- < C x it suffices to restrict 



A 

the sup- norms to the compact interval [0, ( 6+3g *) ]. Since continuous functions on 

compact intervals are uniformly continuous, it suffices that for fixed t the map /C A 9 
f°° dp 

f / j G R+ is continuous. Since the integrand is 

bounded by (^) 2 , we can for A > restrict the integral to a compact interval 
[0, A 2 (e)]. Since the Hilbert transform is continuous according to the discussion in 
( |5.31|) and rational functions of continuous functions with denominator separated 
from zero are continuous, the whole function fC x 3 f i— > -jr= G 

f(p) J 

R+ is continuous for every p and hence uniformly continuous on [0,A 2 (e)]. This 
implies \\T\f — 7a/ < e for /, / sufficiently close in /C A . 

(K2) and (K3) are the conditions in Proposition |5.13| which imply that T A (/C A ) C Cq(R + ) 
is relatively compact. Hence T A (/C A ) C Cq(R + ) is compact, with T A (/C A ) C /C A by (Kl) and 
(K4). This means that T A maps /C A into a compact subset of /C A . The properties (K5) and 
(K6) verify the other conditions in the Schauder fixed point theorem (Proposition [5.11 
which thus guarantees that the map T A : /C A — > /C A has a fixed point Gbo'- 

Theorem 5.15 For any A > 0, the equation 



1 + 6 r \ 64tt 2 



«/ 



has a solution G^o G Cq(R+). This solution is automatically smooth, monotonously de- 
creasing with uniformly bounded derivative — (j^ + C^G^o < < and pointwise 
bounded by < G b0 < j^. At b = one has G 00 = 1 and ^\ b=0 = -(1 + □ 



In retrospect the Theorem justifies the assumption ( |5.9[ ) on Holder continuity of Gab- 
It is clear that by restricting to [0,A 2 ] and ignoring the behaviour at oo, also (|5.27|) 
has a solution G^ G C 1 [0, A 2 ]. Note however that this proof works for fixed A 2 (including 
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A 2 = oo) but without control on the limit. This means that our existence proof does not 
imply the (highly plausible) existence of the limit liniA^oo G^o an d the equality with G^q. 



The solution Gw of ( |5.38| ) provides $b(a) via ( |5.16b| ) and then gives the complete 



two-point function G a b according to ( |5.19|) . Moreover, as we show in the remaining part 
of the paper, all higher correlation functions are explicitly given in terms of Gbo- 

If we also knew uniqueness of the solution, then the resulting unique $0 would be the 
only candidate for a solution of ( p.26| ). But before addressing the uniqueness question 
one has to understand the possible non-trivial solutions of the homogeneous Carleman 



equation which we ignored by Assumption |T2|. This is a project of its own 



6 The planar regular sector 

The remaining part of the paper is devoted to higher correlation functions. In this section 
we prove an algebraic recursion formula for all functions with B = 1. Functions with 
B = 2 punctures are treated in Appendix ||, where also an outlook to B > 2 will be 



given. Subsection 3.3 gives the proof that the /3-function of the model vanishes non- 



perturbatively. Subsection 5.4 shows that all non-planar functions are suppressed. 



6. 1 Derivation of the Carleman equation 

Specifying (|3.9| ) to the parameters and index sets ( f4.3b| ) of self-dual noncommutative 
4 -theory, the planar regular (N > 4)-point function satisfies the equation 

^(0) ^(0) r (o) r {o) 

lr>(°) \ ST \pbi-bN-i\ ^\pbi\ lafei— feiv-il 



z ^lii.-.&jv-ii ~ A Yl 



!(bl - \a\) 

^ G {0) - G (0) 

\r<(°) r>(°) \hi+i—kN-ib2i\ |fe2i+i---fcjv-isl ir~\\ 

As before, the Appoint function G^ bi bN , only depends on the 1-norms \a\ and \bk\ of 

its indices, which allows us to substitute X] P eN 2 f(\p\) ^ S|d=o(W + 1)/(H)- ^ n or der 
to prepare for the limit (|5.1|) and the passage to the variables (|5.2| ) we redefine 



G abl .., N ^ ■= ^^e^-'G^..,^ , (6.2) 



which for A^ = 2 is compatible with ( |4.7|) and ( |5.5|) . Multiplication of (34.) by 
then leads in terms of 05.21) and (16.21) to the formula 



(27re)V 2 G (0) 



Z- 1 A / 4 \ 2 \^(\ 1 u G P^- b ^ ~ ty^-^-i 

G a b! bl '" bN ~ 1 f^W 2 ' ) ^(H-N) 

N-2 

2 /^i s~~1 
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where p,a,bk are viewed as functions of \p\, \a\, \b_k\ according to (|5.2|). In the limit ( |5.1| ) 
this equation converges to an integral equation which under the assumption that G a b 1 ...b N _ 1 
is Holder-continuous can be rearranged as 

Grabi...b]y_i ( Z 



G 



abi 



A 



V-2 
2 



j 

o 2 ; - a 



(6.4) 



Using the previous definitions and identities ( |5.12b| ), fl5.12a| ) and ( |5.7| ), the prefactor 
involving Z~ l is treated as follows: 



1 + 6i^a[C» 



a aG„ 



+ &iG 



A 



oO 



1 + y 64vr 



aG, 



h , i + ^[g.o] 



aG 



iTy- 1 + 6S M »[ ( -- a)G -°] ' (6 ' 5 > 



=0 



Inserted back into ( |6.4|) we obtain for (aG a b 1 ...b N _ 1 ) again a Carleman equation 

A 



— + 

a 



aG, 



64vr 



%a[* G, bl ... bN _ 1 ] 



A 



JV-2 
2 



1 + yy E ^...ba 



Gb 2l+1 ...b N _ 1 b 2l Gb 2l+1 ...b N _ 1 g 

hi - a 



(6.6) 



Remarkably, the Carleman equation ( |6.6| ) has the same homogeneous part as the equation 
( |5.12a| ) for the two-point function D a b 1 , only the inhomogeneity is different. Its solution 
is given by Proposition RTTl again under the assumption G = 0, as: 



JV-2 
2 



G 



ab\...b N -\ 



64tt ^ 



sint?b 1 (a) /Gfe ^, JV _ 1 f ) G& ; +1 ...6 JV _ ia , , 
^, ( ; cos($ 6l (a)) 



e -^[<V] s i n (tf 6i (.)) 



&2/ - a 

Gfe 2(+1 ...6 JV _ 1 (, 2i — G6 2!+1 ...6 JV _ 1 , 



>2Z 



(6.7) 



It gives iteratively the solution for any planar regular correlation function G a b 1 ...b N _ 1 - 



6.2 An algebraic recursion formula for planar regular N -point functions 

The generating functional W[J] for connected iV-point functions is real. Since both 
Gb ...b N -i is a real function of bi G [0,A 2 ] as solution of the real equation ( |6.6|) and the 
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source is real, i.e. J = J* or J pq = J qp , this implies for the original regular iV-point 
function 

^ ] Gp ---PN-l^POPl ' ' ' JpN-2PN-lJpN-lP0 

P0,--;PN~l^jf 

Gp ...ptf^ 1 Jp PN-lJpN-lPN-2 ' ' ' JpiPO ■ 

Po,—,Pn-i£®%- 

Therefore, the regular iV-point function is not only invariant under cyclic permutations 
but also invariant under reversal of the order of its indices (or orientation reversal), 

Gb b 1 ...b N _ 1 = G ? ;, b JV _ 1 ...;, 1 . (6.8) 
Reversing the order of the indices b±, . . . , fojv-i in the Carleman equation (|6.6| ) leads to 

l~ + a~G~ n ) y aG « b »-^) ~ ei^« L G -^- 1 -^J 



JV-2 
2 



A Gi )2l _ 1 ...b 1 b 2 i G\ )2l _ 1 ...b 1 a ,„ „s 

If we now take the difference between ( 6.6| ) and ( |6.9D , all terms which involve the Hilbert 
transform cancel due to the reversal symmetry ( |6.8|) , so that after division by &i — 
and reversal and cyclic index permutations we arrive at 

Theorem 6.1 The planar regular connected (N > 4) -point function of self-dual noncom- 
mutative <p\-theory is obtained from the algebraic recursion formula 



JV-2 
2 



p _ ( ^) Gb Q b 1 ...b 2 i- 1 Gb 2l b 2l+1 ...b N - 1 Gb 2l b 1 ...b 2 i-iGb b 2 i + i...b N - 1 ,„ „x 

6obl - bAr - 1 " (i + y? ^ - b ^ - 6 ^-0 ' 

where 1 + y = — ^jjf~\b-o' The result holds independently of Assumption \5. <\ . □ 



One may speculate that for large indices a,6> 1 there is a scaling relation G sajS b = s r 'G a b- 

■ - N — 4 

Such a relation would by Theorem |6JJ result in G S 6 ,..,s6 jV _i = s^~ v Gb ...bf, 



ojv-i • 



For N e {4, 6} the general formula (|6.10|) specifies to 



^ ( — A) G a bG c d — GadGbc ((\-\-\\ 

abcd ~ (1 + y) 2 (a-c)(b-d) ' 1 j 

p (—A) [ Gb^Gb^bjbs — Gb^Gbzbjbsbo Gbpbjb^sGb^ — Gbjb^bjGb^o 

b0 - b5 ~ (i + yy\ (bo-hXh-h) + ( &0 _ fo4 )( &1 _ fo5 ) 

( — A) 2 f Gb b!Gb 2 b 5 Gb 4 b 3 GbobsGb^Gb^bs , Gb^Gb^Gb^bx 



^02^24 ' ^15^53 ^40^02 ' ^53^31 ^24^40 ' ^31^15 

^ ^ ^ / 1 1 1 1 1 1 

+ ^6061^263 ( 7 — r — ' 7 — 7 1" 7 — 1 — ' 7 — 1 1" 



^02^24 ^13^35 &04&42 ^13^35 ^04^42 ^15^53 

+ G bob5 G b2bl G Hb3 [j—r- ■ T—f— + j—r— ■ j—r— + j—r- ■ 7—7— ) }" > (6-12) 

\O04O42 O53O31 O02O24 O53O31 O02O24 O51O13, 
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where bij 



hi — by The denominators arising in fl6.12| ) are not unique. The identity 
1 1 1 



/ j + nr + nr = ( 6 - 13 ) 

DijOjk OjkOki OkiOij 

can be used to rearrange the denominators. In Appendix |A| we give an alternative proof 



of ( |6.11| ) by direct solution of ( |6.7| ) using Tricomi's identities ( |5.15b| ) and ( |5.15a|) and the 
explicit formula ( |5.19|) for the two-point function. 

We introduce a graphical representation for a contribution to Gb ...b N _ 1 - We put the 
indices b , . . . , b N _i in cyclic order on the circle. A factor G^b is drawn as a thick chord 
between bi and and a factor , 1 , as a thin arrow from b m to b n . Reversing the arrow 
thus changes the sign. Later on we connect the arrows to trees where it will be necessary 
to specify the root as a thick point. The corresponding graphical description of ( |6.11p 
reads 



G 




(6.14) 



The two graphs are turned into each other by a rotation of =j- and possibly reversion of 
both arrows. 

Our rules give the following graphical representation of the six-point function (|6.12 ): 



G h 



two rotations 

2-Kk 

k e {1,2} 





one 
+ rotation 

2tt 



y (6.i5) 



by 



6 



We first notice that only the 5 non-crossing pairings of the cyclic indices bob^b^b^b^ 
appear. The remaining pairing Gb Q b 3 Gb 1 b 4 Gb 2 b 5 has a self-crossing. The non-crossing 
pairings of 2n cyclic indices are counted by the Catalan number C n = 



(2n)! 



n!(n+l)!' 

Note that the 4 possible products of even/odd paths starting with root b /bi are 
distributed over the pairings in such a manner that the paths intersect the pairings only 
in the vertices and every product of paths appears exactly once. Further rotation of the 
first line by pp and of the second line by pp give the same graphs after common reversal 
of the arrows and use of the identity (|6.13| ). 

These observations generalise to any planar regular iV-point function: 



is a sum 



with pref actors ± ( ) 



2 " 



Proposition 6.2 The N -point function Gb ...b N _ , -••<■.....,,, _■_ v ( { ^ y)2 

of graphs with vertices bo, . . . , fr/v-i put in cyclic order on a circle and (fat) edges between 
bk, bi representing a factor Gb^ and (thin) arrows from b^ to bi representing a factor b _ b . 
We call b 2 i an even vertex and 6 2 «+i an odd vertex. In every such graph we have (or can 
achieve via (\6.1ty): 
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1. Every even vertex is paired (by a fat chord) with exactly one odd vertex, and these 
chords do not cross. 

2. Taking any even vertex as a root, there is a an oriented rooted tree between all even 
vertices. This even tree has no self-intersection and intersects the fat chord only at 
the vertices. 

3. Taking any odd vertex as a root, there is a an oriented rooted tree between all odd 
vertices. This odd tree has no self-intersection and intersects the fat chords only at 
the vertices. 

Proof. By ( |6.10| ), the iV-point function Gft .- &iv-i ar i ses by connections of 2/-point func- 
tions with (N — 2/)-point functions via an even arrow and an odd arrow. We proceed 
by induction and assume validity of the Proposition for the 21- and (N — 2/)-point func- 
tions, starting with Gb bx which is represented by one chord without any arrows. We can 
symbolise ( |6.10| ) (without prefactor ) as follows: 




Here, the * locates the vertices and the dashed line separates the 2Z-point function from 
the (N — 2/)-point function of which we do not show their inner chords and trees. 

1) By induction in one of these smaller functions every even vertex is paired with 
exactly one odd vertex. These pairings do not cross in the smaller function and 
since they do not cross the dashed cut, there cannot be a crossing between chords 
of the two smaller functions. 

2) All even vertices in the 2Z-point function and in the (N — 2/)-point functions are 
connected by a tree which by induction has no self-intersection and intersects the 
chords only in the vertices. The additional arrow from 6 to bi connects both trees 
of I and — — I vertices to a larger tree of y vertices. It is geometrically clear that the 
additional arrow intersects both the two smaller even trees and the chords precisely 
in the vertices b and b\. 

3) All odd vertices in the 2Z-point function and in the (N — 2/)-point functions are 
connected by a tree which by induction has no self-intersection and intersects the 
chords only in the vertices. The additional arrow from b\ to bx-i connects both 
trees of I and y — I vertices to a larger tree of y vertices. This larger tree does not 
have self- intersections, but in general intersects the chord emanating from bo- Since 
a cyclic permutation of the indices interchanges even and odd vertices and hence 
even and odd trees, by 2) it must be possible to use the identity ( p,13| ) to rearrange 



the odd tree in such a way that the possible intersection with the chord emanating 
from bo can be avoided. □ 



39 



Proposition p.2| determines the planar connected iV-point function up to the problem 
to identify the product of even and odd trees associated with a given non-crossing chord 
diagram. Due to (|6.13|) this product cannot be unique. We leave this question of canonical 
trees for future investigation. We finish this subsection with the eight-point function which 
shows that not all possible trees do actually arise: 



Gb Q ...b 7 

( ^) ( Crf, b 1 (j6 2 63 64 65 65 67 G b\b 2 G 636465656760 f G b bib 2 b 3 G b 4 bc,bgb 7 G b^b 2 b 3 b 4 G b^b^b-jb^ 



G 6q6i 62636465 C 6567 ^616263646566^6761 



+ 



(60 -64) (6l 



(&0-&6) (6l 



:-a) 



Gb b 1 Gb 2 b 7 Gb 4 b 5 Gb a b 3 ( 



+ (j6 6 3 G r 6 2 6 1 (j6 4 6 7 (j6 6 6 5 {b 20 b Q4 b i6 

+ Gb b 5 Gb 2 b 3 Gb 4 b 1 Gb 6 b 7 {b 24 b 40 b 06 
+ Gb Q b 7 Gb 2 b b Gb 4 b z Gb 6 bi ( 642626^60 

+ Gb biGb 2 b 7 Gb 4 b A Gb 6 b 5 {b 02 b 24 b 4 e 

+ Gb b 3 Gb 2 b 1 Gb 4 b s Gb 6 b 7 ( 620604646 

+ Gbob-jGb^Gb^Gbebs ( b 24 b 46 b G0 

+ Gb biGb 2 b 5 Gb 4 b 3 Gb 6 b 7 {b 42 b 26 b ao 

+ Gbab-rGb^Gb^Gb^b-i {b i6 b fi0 bo 2 

+ Gb^Gb^Gb^Gb^ {b 64 b 40 b 02 

+ Gb b 5 Gb 2 b 1 Gb i b3 Gb 6 b 7 ib 60 b Q2 b 24 

+ Gb b 7 Gb 2 b 3 Gb 4 b 5 Gb 6 bi {b oa b 62 b 2 4, 

+ Gb^Gb^Gb^Gbeb, { bo6 bl 4 b i2 

+ 1 1 + 

602624646 613635657 

+ 1 1 + 



602626654 617673635, 

1 ) 

613637675 / 

1 ) 

631615657/ 

1 ) 

635651617 / 



1 
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1 


1 


+ 


1 


1 


617673 635 


602624646 


617675 653 


602626664 


617675653 


1 


+ 


1 


1 


+ 


1 


1 


613635657 


620606654 


613635657 


620606664 


613637675 


1 


+ 


1 


1 


+ 


1 


1 


631615657 


624646 660 


631617675 


624640606 


631617675 
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1 


1 


+ 


1 
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635657671 


642620606 


635657671 


642 620 606 


635651617 


1 


+ 


1 


1 


+ 


1 


1 


653637671 


645660602 


653631617 


646 662620 


653631617 


1 


+ 


1 


1 


+ 


1 


1 


657671613 


664642 620 


657671613 


654642620 


657673631 


1 


+ 


1 


1 


+ 


1 


1 


675651613 


660 602 624 


675653631 


660 604642 


675653631 


1 


+ 


1 


1 


+ 


1 


1 


671613635 


605664642 


671613635 


606654642 


671615653 


1 


+ 


1 


1 


+ 


1 


1 


617675653 


606654642 


617673 635 


605654642 


613637675 



602625654 613635657 

1 1 



+ 
+ 



606662624 613635657 

1 1 



606662624 617673635 606662624 

I si r< r< ( l l I l l 

^^6067^6261^6463^6565^606664642 675653631^606652624 675653631 

1 1 + I 1 + 1 



613637675 
+ 1 



+ 



1 



606 604642 
1 



613615657 / 



602625654 675653631 
1 



+ 



602624645 671613635 
1 1 



+ 



602624646 671615653 
1 1 



602626664 671615653 606662624 675651613 602626664 675651613 



+ 



602624646 675651613 
1 1 



+ 



602604646 675673631, 

(6.16) 



We have used (|6.13|) to bring the denominators for the last two products of Gb ibj into a 
form which (in the sum) is symmetric in the even and odd indices and has a common 
sign. A novel feature compared with the six-point function is the appearance of trees with 



branches 



and 



, , . Trees become unavoidable in the 12-point 

5673631 r 



606604642 613615657 602604646 6756 
function where some pairings do not permit a non-crossing rooted path. It seems natural 

to require that paired roots carry the same number of branches. If we furthermore require 

that in 8-point functions only the root has more than one branch, then the representation 
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given in ( |6.16p seems distinguished. Its graphical representation reads 




Again we obtain all C4 = 14 non-crossing pairings of the cyclic external indices b , . . . ,b 7 
which give the products of G^. . 



6.3 The effective coupling constant 



Of particular interest is the effective coupling constant X e jj = — Goooo- Indeed, we have 
\ e g = r ooo where the 1PI function T abcd is obtained by amputation of the connected 
two-point functions from G abc d and a change of sign, i.e. T abcd = — „ f? ah A d „ ■ The 



GabGbcG c dGda 

direct computation from (|6.11|) involves difficult limiting processes. We therefore compute 



Goooo = lim a ^ G a00 oo directly from flfTTD , using G 00 
cos ^o(O) = sign(A): 



1, lim 



— lim 



647rsini9o(fl) (G 



lim 

a->0 



a->o (1 + y) 2 a 

647r sin $o(a 



■cos(0 o (a)) + e<^H A a 



O-5-0 



-n A .m 



GgQ— 1 



sin(# («)) 



-(l + y) and 
G.o-1 



[i + y)a 



sign(A)+7^ 



o H A lM~H A .lM 



sin(# («)) 



1 — G.o 



Recalling lim 



a->0 



sini?o( a ) 



64?r 



from ( |5.25| ) and Lemma |5.3| , we obtain 



A 



64vr 



y 



:^+y) 2 

41 



n/A 
tin 



1-G.o sin 2 tf ( 



G,o 



Spelling out the Hilbert transform, taking ( p.2U| ) into account and going to the limit 
A — )• oo, we have proved: 



Proposition 6.3 The effective coupling constant \ e ff = — Goooo of self-dual noncommut- 
ative <fi\-theory is given in terms of the bare coupling constant A by the following equivalent 
formulae 

> A 64A 2 f°° l-G p0 sm 2 Mp) 

X eff = + M , ^ / dp 



l + y (l + ^) 2 7o PG P0 (Ap) 2 

1 — G p o 



Theorem \5 . 1 b\ the change X e jj ^ X is only a finite renormalisation of the bare coupling 



constant in response to an infinite change of scales, which means that the QFT model has 
a non-perturbatively vanishing (5 -function. □ 

We recall that vanishing of the /3-function at any order of perturbation theory was 
proved by Disertori, Gurau, Magnen and Rivasseau in |pGMR07|| . A decisive step in 
their proof was to neglect a term in |DGMR(T7| , eq. (4.18)] which in our notation reads 



^p eN AA W\ ^pmOm' This is perturbatively justified, but as far as we can see, this argument 
cannot be used non-perturbatively; at least it is not obvious. Our derivation of (|6.18|) 
shows that the omission of the mentioned term is not necessary in the non-perturbative 
treatment. Of course we have used the fact that the /3- function vanishes in our ansatz for 
the action functional. If instead of (|6.18| ) we had arrived at a relation in which Goooo di- 
verges, then indeed it would have been necessary to replace the coupling constant A \Z\ 
and to eliminate the ill-defined Z\ by the normalisation of Goooo- This procedure might 
be necessary for other models. For self-dual noncommutative 4 -theory in 4 dimensions 
we do not need to renormalise the coupling constant. 



6.4 Suppression of the non-planar sector for 6 — > 00. Ward identity revisited 

Inserting the scalings ( p.2[ ) of the one-cycle functions and ( |B.8[ ) of the two-cycle functions 
into we see that some parts of CrKi are suppressed over G$ b , in the thermodynamic 
limit 9 — > 00. We prove that this holds in general: the whole non-planar sector with genus 
g > 1 is scaled away for 9 — > 00. 

Proposition 6.4 The following rescaling 

g^m*.*. = ^(^ 2 > 2W,+ +Ni, - 4+2fl+ Nl i v 1 i..«...>y • <« 9 > 

of the (Nx+ . . . +Nb )-point function of genus g has a finite limit for 9 — > 00. 

Proof. We proceed by induction, using the fact that the planar sector g = is already 
verified in previous sections and Appendix [FJ. We start with the Schwinger-Dyson equa- 
tion ( |3.6| ) for G a b and use E a and A 4 given by ( |4.3b| ). Multiplying ( p.6| ) at g > 1 by 
fi~ 2 (27r9jj 2 ) 2+49 , and using that /1 2 E and ^^jFf 2 nave limits, it is straightforward 
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to check that both sides of ( |3.6| ) have the same scaling ( |6.19| ). The same considerations 
prove ( |6.19|) for the iV-point generalisation ( |3.7| ). The only difference is the last term 



G 



(a 1 ) 

\bi...b 2l \' 



G 



(s") 



G 



(a") 



>2l + l'-- b N~l b 2l\ \ b 2l + l--- b N-l a 

Si,„, — E a 



which is readily checked. 



These building blocks reappear in the Schwinger-Dyson equations for B > 1. By 
universality of the scaling behaviour (see also | GW05b |) one confirms ( 6.1 9|) in the general 
case. □ 



^\b\...b\\...\bf...b* B Y 



Recall that the genus expansion reads G[ 6 i_ 6 ^ |...|&b...&b | = 

Expressing G™' in terms of the functions G 9 on the lhs of ( |6.19D which have a limit for 
9 — t- oo we see that a genus-g contribution is suppressed by a factor (27i9fi 2 )~ 4 ' 9 . Hence, 
the limit 9 — > oo suppresses the entire non-planar sector with g > 0. Note that ( |6.19| ) also 
suggests that iV-point functions with B > 1 are suppressed by a factor (27i9fi 2 )~ 2( - B ~ 1 ^ over 
the one-cycle function. However, planar iV-point functions with different B contribute 
to the logarithm of the partition function (in contrast to the genus expansion) with a 
different structure of external sources. Therefore, we are not allowed to compare these 
functions. We recall that the 9 — > oo limit of the model at O = was studied by Becchi, 
Giusto and Imbimbo |BGI02| , |BGI03|| . They already established the suppression of the 
non-planar sector (also noticed in | MVS00 |) and presence of functions B > 1 ("swiss 
cheese") in the limit 9 — > oo. 



Functional derivatives of the original Ward identity 



express the index integral of 



an iV-point function in terms of other functions. Putting p = &2j 



6jv and applying 



the derivatives with respect to 



we obtain from (EJB 



{E b2 - E bN ) ( ^ G\nb 2 ...b N \ + ^2 G \b2-b k \b k ...b N \ + ^ G \ b 2- b 2l+i\ G \b2l+i-b N \ 



nel k=2 
G\b N b 3 ... bjv_i| — ^|6 2 63— bjv-i| • 



(6.20) 



Inserting 3 bp and the rescaling ( |6.19| ) we see that the functions with B = 2 cycles 
increase the genus by 1. The restriction to the planar sector g = thus leads with (|5.2|) 
to the equation 



z(i + y)(b 2 -b N )((i + y) 



pdp 



N-2 
2 



G p b 2 ...b N + ^ ] Gb 2 ...b2i +1 Gb 2 i +1 ...b N 



1=1 



/ 64vr 2 

= G riN b 3 ...b N _ 1 — Gb 2 b 3 ...b N _ 1 ■ (6.21) 
It is interesting to verify (|6.21| ) for N = A. From the explicit formula (|6.11|) we obtain 



Z ^^ X (G b cU A c [*G. d ) - G dc H K c [*G.b}) + Z(l + y)(b - d)G hc G cd = G dc - G bc , (6.22) 



in agreement with 
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7 Conclusion and outlook 



In this paper we were able to non-perturbatively solve a toy model of four- dimensional 
quantum field theory Our solutions of the correlation functions can be viewed as summa- 
tion of infinitely many renormalised Feynman graphs (see also Appendix The main 
tools for these achievements are Ward identities in field-theoretical matrix models, the 
resulting Schwinger- Dyson equations and the solution theory of the Carleman singular in- 
tegral equation. The 2-point function plays a central role in our approach. It is expressed 
in terms of a single function G of one variable given as solution (that we proved to exist) 
of a non-linear integral equation. We have found remarkably simple algebraic formulae 
for the higher correlation functions which are expressed as weighted sums of products of 
2-point functions. 

We are left with a few problems: 

• We have ignored so far the possible non-trivial solutions of the homogeneous Car- 
leman equation. The numerical treatment in [ pW13|| justifies this for coupling 
constants < A < 6Att. This needs a rigorous confirmation, which is only possible 
if the whole space of solutions of the Carleman equation is taken into account. The 
investigation of the consistency requirements should then select the right point in 
the solution space. In this way it should be possible to continue the solution of the 
model to A > 64tt. It seems plausible that the true solution changes its sign for 
A > 647T. This could be interpreted as the impossibility of an infinite correlation 
length so that A = 64-7T might be the end point of the A-family of critical models. 

It makes little sense to prove the (important!) uniqueness of the solution of the 
integral equation for the two-point function before clarifying the freedom from the 
homogeneous Carleman equation. 

• Clarify the combinatorics of the weight factors in the formulae expressing the Ap- 
point function in terms of the 2-point function. We already know that the non- 
crossing chord diagrams arise which are counted by the Catalan numbers. Identify- 
ing the combinatorics might suggest links to other models. We have already checked 
that the same combinatorics arises in the 2-dimensional variant of this model; only 
the 2-point function is completely different. We have also some evidence that higher 
correlation functions of the original discrete matrix model, i.e. for finite 9, satisfy 
the same combinatorial formula. 

• Reformulate the infinite volume limit 9 — > oo in position space. On a formal level, 
this limit restores translation invariance, which needs a rigorous confirmation. It 
would be interesting to know the resulting interaction term. 

• If the position space reformulation succeeds it seems not impossible that the model 
satisfies reflection positivity. Although it cannot be a model of usual relativistic 
quantum field theory (because locality is badly violated), a certain variant could 
result which deserves further investigation. 

• We feel that getting thus far with this model was only possible because there is a 
deep mathematical structure behind. It could be a sort of integrability or struc- 
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tures similar to those in conformal field theory. Searching for this structure seems 
promising. 



A Another proof of the formula for the four-point function 



It seems surprising at first sight that we were able to compute the iV-point functions 
purely algebraically, i.e. without the analytic solution of the Carleman equation (|6.6| ). To 
remove any doubt with this procedure, we supply here the proof that the direct evaluation 
of the Carleman solution ( |6.7| ) leads after insertion of the explicit formula ( |5.19| ) for the 
two-point to exactly the same result ( |6.11D for the four-point function. 



Lemma A.l The following Carleman equation 

Aa n/A 



a aG a0 



abed 



^ T? 1 - ^cd G ad 



64tt 



c — a 



has the explicit solution 



abed 



(a — c)(b — d) 



G 



GabGcd 



ad 



G 



be 



(A.la) 



(A.lb) 



Comparison of ( p.6| ) with ( |A.laj ) shows aG abcd = Gb c Fabcd- Therefore, we reestablish 

the algebraic formula ( |6.11| ) for the planar regular renormalised four-point function G abod . 

Proof of the Lemma. The Carleman formula flSJjj ) and rational fraction expansion gives 
the solution as 



abed 



64vr . 

— — sin v b {a) 

A 

647T sini^a) 
A c — a 



647T sini^fa) 



G 



c,d 



c — a 



c 



sm(M*)){Ged-G. d ) 



{Ged-G ad ) cos($ b (a))+e^H A a 

- e^H A [e~*M frinOW) - G. d ) 
G cd cos(Mc))e nAM - n ^ b] - G ad cos(Ma)) 

-h^m siri (^(.))G., 



A c — a 

- e^H A sin(tf 6 (.))G J + e^H A c 



(A.2) 



The last identity follows from ( |5.15a ). The terms H A 



d sm{$ b {.))G. d 



in the last 



line are computed by insertion of Theorem [5?7| for G. d and rearrangement of the resulting 
product sin$<z(») sin$6(«) by (|5.18|) : 



64vr 



e -«fc[*o]^ 



-H A . sin(i?fc(«)) sin(fl d («)) 



w 

sign(A) n A [M A 
b-d 



,H A .[# d -d b ] 



sin($ 



#*{•)) 



sign (A) 
b-d ' 
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The last identity is ( p,15b| ). Inserted back into ( |A.2| ) we get 

F abcd = ™* A SinMa) (G c ,cos(^(c))e^^-^W - G ad cos(Ma)) 
A(l + y) c — a \ 



b-d 

The exponentials of are 



U Uj 

rational functions of G 



si 



, sig: v 
+ — -e 



now expressed by Theorem |5.7] as 



and smd: 

647Tsin^ fe (a)/ aG a6 sin $ fe ( c ) 

^afecd = — : G cd cos(i7&(c)) — — ; — — - - G ad cos($ b (a)) 

A c — a V cGcftSini/f^a) 



AaG a d 



a 



- cos 



(^(a) - 6 (a)) 



64vr(6- d)sin0 dv ~< 

AaG afe G cd sin$ fe ( c ) , ,\ 

647r(o — a) sm 175(a) G c 6smi>d(cJ v ' J 

1 \G ah G cd ( 647ia(b-d) . 

= 7 \7Z ~nl — 7^ \ cos(^ fe (c)) sin i? 6 (c 

(a — c)(b — d) I G 6c V A c 

sin^ fe (c) , s\ 

smiy d (c) / 



/ 647T sin $fo (cl) s 

■ G ad (^—(b - d) cos($ b (a)) sin(i? b (a)) + a- — —— cos {■d d {a) - $ b (a)) 
V A smv d \a) 

It remains to insert (|5.18 ) in order to eliminate the terms ^P-(b — d) in favour of ratios 

Then standard trigonometric addition theorems yield the 

□ 



of trigonometric functions 
assertion. 



B Correlation functions with two boundary components 

B. 1 Two cycles of odd length: Schwinger-Dyson equations 

The case of B = 2 cycles needs distinction of several cases. For the (l+l)-point function 



G\ a \c\ 
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(E a + E a )Z[0] 



^apfpufna 



d{-V) 
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Z[J] 
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(-A 4 ) cd*{{Wl[J} + WZ[J])Z[J]) 
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dJaaddcc 



p,n£l 



E p — E a dJ pa dJ t 



d 2 ( <9Z 



pa^^cc N OJ an OJ n p/ J J=0 



-A 4 ) 



G| a |c| ^C|a|a| + ^ G| 
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+ G| a | a | a | c | + ^G| a | c | an | + G| c | aaa | + G| a | cac | + 2G| a | a |G| a | c | 



46 



| G\g\ c \ - G\p\ c \ | G| ac | - G\cc\ i 
~l E p — E a E c — E a J 

Next we generalise ( |B.1| ) to (B = 2)-cycle (N > 4)-point functions where one of the 
cycles has length 1: 

q tM ( d N ({wM + w 2 a [J])z\J\) 

|a|ci - CAr - 1 ' (E a + E a )Z[0] I dJ aa dJ ClC2 . . . dJ^c^dJc^ 

x - 1 / dZ__ dZ \\ 

p n gj P a P& C].c 2 ' ' ' c N-2£N-i cn-iCi u ^an UlJ np' J J—<J 



The same discussion as for the derivation of fl3.5| ) applies: 

G ! | a | Cl ... Cjv _ 1 | = — | ^C|a|a| + ^ ] C|an| ^ G ? |a|ci...c JV _ 1 1 ~t~ ^|a|a|a|cj....cjv-i| + G\ ac ,, , , v , 



iV-1 

+ G ? |a|an|ci..,Cjyr_i| + C| a | Cl ... CfeaCfe ... Cjv _ 1 | + 2G| a | a |G| a | Cl ... Cjv _ 1 | 

nel k=l 



E 



G ? |p|ei...Cj V _i| C| a | Cl ... Cjv _ 1 | 



-En — E, 

N-l 



pel - p " a 



EG\ Cl ...c k -ic k c k c k+1 ...c N _ 1 \ G\ Cl ...c k _ 1 c k ac k+1 ...c N - 1 \ "1 

J ' ( ) 

fc=l Cfe a 



It is not difficult to generalise this formula to the case of two cycles which both have odd 
length. To simplify the notation we write 777 j£- — ttj =: » lk k 



G\ab 1 ...b 2 i\c 1 ...c N _2i- 1 \ 



>PlP2--- J P k -lP k uu P k Pl vu Pl—Pk 

:-A 4 ) f d N {(W^[J} + W![J})Z[J]) 



(E a + E b )Z[0]i dJ^dJt^ 
8 



d N iJ — -J! 9Z 

^ I " T)T). ^ -r tJ 1 



dJ an dJ n p 



^ ^ E p - E a dJ pbl dJ blb2 ■ ■ ■ dJ h2l _ xh2l dJ h2ia dJ^-} 2i i ) j=o 



The derivatives are evaluated according to the discussion preceding (|3.5|) . We give this 
formula directly in its genus expansion, obtained by similar considerations as in Sec. 

\ab 1 ...b 2 i\c 1 ...c N _2i- 1 \ 

_ (-A4) f \^ n {g>) n (g") . \- \- n {g') n {g") 



f G (9') G (g") , G ( 

(_ / j |afei...6 2 il c i--- c JV-2i-il \a\a\ ' / / / -> \ab l ...b 2 i\c 1 ...c N _ 2 i- 1 \ [< 



ft _|_ ft \ ' / |afei...b 2 il c i--- c JV-2(-il |a[a| ' / / / / |abi...fe 2 il c i--- c iv-2i-il |<H 

a 1 9'+g"+l= 9 ne/g'+ 9 "= 9 

|a|a|afei...6 2 il c i--- c iv-2i-il |aaabi...6 2! |ci... cjy_ 2 j-i| / > |an|a6i...6 2 i]ci---cjv— 2Z-i] 

2« N-21-1 

Z— ✓ |6i...bA.o6 )fe ...62(a|ci...c iV _ 2i _ 1 | ~ r Z-e |ci...c fc oc fc ...cjv_ 21 _ 1 |o6 1 ...6 21 | 
fc=l fe=l 
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q(9) _ q{9) 

' 2-^i |a|a&i...& 2 j| |o|ci...cjv_2t_i| 2-^i E — E 

g'+g"+i=g pel p a 

1 6i ...6 fc 1 6 fc _,_i ...6 2 z bfc I ci ...cjv— 2i— 1 1 l&i— h\h+i— 62ja|ci...cjy_2i_i| 



E 



Eh, — E„ 

k=i bk a 

i r (g") _ r (g") 

E\ ^ Q(g') ^\b 2j ...b 2 lb 2j -i\ ^\b 2j ...b 2l a\ 

|6i...6 2 j-i|ci...c J v_2(-i| TP _ TP 

g'+g"=g i=i " 2j - 1 

G (g") _ G ig") 

(g') |b2j+i— fc2ifr2j-i|ci— Cjy_2i-i| |6 2 i+i...6 2 ia|ci...cjv-2(-i| 

■ &2 J'I Ek —E 

I i // '1 »2i a 

g'+g"=g j=i 3 

AT-2Z-1 _ £f(s) 



E E«E 



E"|c 1 ...c fc 6i...62 ! c fc c fc+1 ...c JV _2!-ll |ci...c fc 6i...6 2i ac fe+1 ...c ]v _ 2 i_i| 1 / D Q \ 
r ' (B3) 



fe=i 



For ! = we have to put E —jjr (-> E ^ E . The restriction to the planar (g = 0) sector 
gives: 

Proposition B.l Given a real ^-matrix model with infective external matrix E. Then 
the (unrenormalised) planar connected N -point functions with two cycles of odd length 
satisfy together with the single-cycle functions (\3.8[) and (\3.Q) the system of equations 



q(0) q(0) _ q(0) q(0) 

(0) _ ^ |pbl---M c l--- c iV-2i-il \abi\ |a6i...b 2i |ci...c JV _2i-i| |pf>l| 



^\ab 1 ...b 2l \c 1 ...c N _ 2 i_ 1 \ ^ 4 ^ ] 



i r (o) _ r (o) 

, x \^ r (0) r {0) ^\b2j...b2ib2j-i\ Ur |b 2 j-b2i"l 

+ M 2^ U \ab 1 \ U \b 1 ...b 2j . 1 \c 1 ...c N _ 2l _ 1 \ E _ E 

j—l b 2j-i a 

I q(0) _ q(0) 

,\ V^/-r(0) sy(0) \b2j+i---b2lb2j\ci—c N -2i-\\ \b 2j+1 ...b 2 ia\c 1 ...c N - 2 i- 1 \ 

+ M Lr |a6 1 [ C %...& 1 y| Eh -E 

j=l ^bti 11 ' a 

N-21-1 q(0) _ q(0) 

-i(0) \c 1 ...c k b 1 ...b 2l c k c k+1 ...c N _ 2i _ 1 \ |ci...c fc 6i...fr 2 iacfc + i— Cjy_ 2 i-i| 



+ A4 Yl G \ab 1 



I E — E 



k=l Ck 

(B.4) 

For I = one has to replace by G ^ 

Proof This follows from the (g = 0)-case of ( |B.3| ) and elimination of J2nei ^[oni v * a (HID- 

I ! 



B.2 Two cycles of even length: Schwinger-Dyson equations 

The situation is different in the case of two cycles which both have even length because 
the action of the two J-cycles on Z decomposes into its action on log Z and the separate 
action of each one cycle on one of the factors log Z in |(log2) 2 . This means 

Gr\ a b 1 ...b 2 i_ 1 \c 1 ...c N _ 2 i\ + G^|o6i...6 2 ;_i | G\c!...c N _ 2 l I 
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-A 4 ) <d N ({Wl[J] + WM)Z\J]) 



(E a + E b )Z[0}\ dJX... b2l _dJ»Z 



21 

CN-21 



d N (j —-J 9Z 
\j \ o rm. ^ ^ 



I " V pn f)T " ILLl r)T J ^ 

\ 1 V (J J an (JJ np / 1 

~ii e p ~ E * dJ pbl dJ blb2 ■ ■ ■ OJb^b^dJb^adJ^ 2 ^^ J J=o ' 



In addition to the discussion of ( |3.5| ) there is now the possibility that the derivatives with 
respect to J 2l bi b and J^~.c l N ^ 2l ac ^ separately on °ne of the functions W*[J] + W% [J] 
and Z\J\. The analogous consideration applies to the last line. In this way all terms 
which constitute G\ Cl ...c N _ 2 i\ times the rhs of equation (|3.5| ) for G r [ a 6i...6 2 i-il are generated. 
Cancelling these terms with the lhs we obtain after genus expansion 

q{9) 

|a6i...6 2 !_i|ci...cjv-2il 



I V G {9 ' ] G {9 " ] + V V G ( 



E + Eh I X J |0"l-' ) 2!-l|ci...cjv_2!l ' 1 / j / 4 ' \ab 1 ...b 2l _ 1 \c 1 ...c N _ 2 i\ \an\ 

bl g'+9"+l=9 g'+g"=gnel 

\a\a\ab 1 ...b 2l _ 1 \c 1 ...c N _ 2l \ ' / j \a\a\a ...c N _ 2l | |a6i...6 2 i-l| 

g'+g"+l=g 

/ |on|o6i...6 2! _ 1 |ci...c J v_2il / / / |m|ci...c iv _2il |o6i-62!-il 

ne/ g'+g"=g n£l 

21-1 

+ G {9 ~ 1) + S >^G {9 ~ 1) 

' \aaab 1 ...b 2l _ 1 \c 1 ...c N _ 2l \ ' / j \b 1 ...b k ab k ...b 2l _ 1 a\c 1 ...c N _ 2 ;| 

fc=l 

N—21 N-21 

Z. ' \c 1 ...c k ac k ...c N _ 2l \ab 1 ...b 2 i_ 1 \ / , / , \c 1 ...c k ac k ...c N _ 2l \ |a&i ...6 2J _ 1 1 

k=l g'+g"=g k = l 

Q(g) _ nig) 

E |p&i...b2i-i|ci-cjy_2i| |a6i-b 2 i-i|ci-CAf-2il 
E v - E a 

pel p a 

21-1 /y(s-i) _ Mg- 1 ) 

\b 1 ...b k \b k+1 ...b 2l _ 1 b k \c 1 ...c N _ 2l \ \b 1 ...b k \b k+1 ...b 2l _ 1 a\c 1 ...c N _ 2l \ 



E b , — E n 

k=i ° k a 

1-1 n(g) _ nig) 

E G (g) \b2j+i---b 2 i-ib 2 j\ |b 2 j + i...fc>2i-ia| 

|6l...6 2 j|ci...c JV _2i| TP, _ ft 

Z-l ^(9) _ ^(9) 

1(9) |&2j+l—f>2i-l&2j|ci— CjV-2i| |b2j+l---b2i-lO-|ci...Cjv_ 2 ;| 



• &2 -' 1 Eu - E 

J ^b 2 j 



N—21 (j(g) _ Q{g) 

" / • (B 5) 

The restriction to the planar ((7 = 0) sector gives: 

Proposition B.2 Given a real ^-matrix model with injective external matrix E. Then 
the (unrenormalised) planar connected N -point functions with two cycles of even length 
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satisfy together with the single- cycle functions i\3.8j and (3~1)) the system of equat- 



ions 



G 



(o) 

|o6i...6 2 i-i|ci...cjv-2i| 

A 4^|a&iHa&i...&2!-ll tj |an| Cl ... Civ _ 2i | A ^ ^...b^. 



N-21 
k=l 



(0) 

\c 1 ...c k ac k ...c N _2i\ 



+ A 4 £ 

pel 
1-1 



G 



(o) 



G (o) _ ^(0) 



G 



(o) 



+ A 4Z^ G |afei| G |6l.. 



[p6i...5 2 t_i|ci...c i sr-a!| [o&lj |afei...fe2i-i|ci---CjV~2il \pbi\ 

E p — E a 

q(0) _ q(0) 

\b2j + i---b 2 i-ib 2 j\ \b 2 j+i-~b 2 l-ia\ 



.6<y[ci...cjv_2t| 



+ A 4 2^ (j |a6 1 | G |bi.. 



•fcil" 



N-21 



q(0) _ q(0) 

\b2 j+i-..b 2 i-ib2j\ci...c N -2i\ \b2j+i...b 2 i-ia\ci...c N _2l\ 

Eb 23 - E a 



+ A 4 G \ab 1 



G 



(0) 



(0) 



(0) |ci...Cfe6i...&2i-iCfeCfc + i...cjv_2i| |ci...c fc 6i...6 2 i_iac fc+ i...cjv_2ll 



fc=l 



(B.6) 



Proof. This follows from the (g = 0)-case of ( |B.5| ) and elimination of Ylnei G \aL\ v * a (HID 



□ 



5. 3 Two cycles of odd length: solution 

r \a\c\ 

= where we have to 

we thus have 



The planar (l + l)-point function G^L satisfies the Schwinger- Dyson equation (p.4j) for 
N = 2 and / = where we have to replace by Cf^i. With indices a, c,p £ 1 = 



/o(0) ^(0) _ ^(0) ^(0) (0) ~(0) 

= A4 2. W^E n + Er _ En ■ 



2 p s 



Again, this equation only depends on the 1-norms ( |5.2| ) of the indices. It follows from the 

r \°k\ 



formulae for Eg and A 4 given in Q4.3b|) and the rescaling (2it9) 2 ^ 2 G^\ = G ao that ( |B.7|) 



has an infinite volume limit 9 — > oo iff we define 



G a]c := (27r^)VG^ | • (B.8) 

____ f nnc laaHc ir\ 

(l+y)G aa 



Multiplication of (p.7|) by zf\+y)G ^us ^ ea( ^ s ^° 



Z- 1 G alc A / 4 Vf/un fr-^ A G cc -G ac 
(l + y)G aa 6A^\e^(l + y)J ^ ] P~a (l + y) 2 c-a ' 

|p|=0 

(B.9) 

where p and \p\ are related by ( p.2| ). Now the thermodynamic limit ( |5.1| ) to continuous 
variables a, c, p £ [0,A 2 ] exists. Under the assumption that G a | c is Holder-continuous, we 



50 



can move -^G pa to the lhs of ( |B.9| ), and taking (|6.5| ) at b\ — a into account we arrive at 
the following Carleman equation for the planar (l+l)-point function: 



a aG a o 



(aG a \ c ) 



^ n/A r p i _ A G cc G a 



(B.10) 



We let $(a) := $ a (a) and stress that the (non-)presence of the subscript at $ will be 
important in the following calculation. The Carleman formula (|5.14|) and rational fraction 
expansion gives the solution of (|B.10|) as 

64vr sintf(a) (G cc -G a 



G 



a\c 



[1 + y) 2 a \ < - a 
647r sin-i^a 



cos(0(a)) + e nA ^n A a e-™ sin(tf (.)) 



G CC G 9C 



[l + y) 2 a(c-a) 



(G cc - G ac ) cos(tf(a)) + e™?^ e-*M sin(tf(.)) (G cc - G. c ) 



— e 



e-^sm($(.))(G cc -G. c ) 



64:% siwd(a)e H ^ 
(l + y) 2 a(c-a) 



G cc cos($(c))e- nA ^ - G ac cos(^(a)) ( 



- n r . 



-^[ e - W - M sin(^(.))G. c ]+^[e-^^sin(^(.))G. c ]) . (B.11) 

The last identity follows from ( 5.15a ). We see no possibility at the moment to further 
simplify this expression. A non-trivial consistency check is the symmetry G a \ c = G c \ a . In 
perturbation theory up to 0(X 2 ) we confirm the symmetry and the agreement with the 
Feynman graph expansion, see Appendix |Cj. 

Next we consider the Schwinger-Dyson equation (|B.4|) for N > 2. It follows from 



QObp and the rescaling ( |6.2|) that ( |B.4|) has an infinite volume limit — > oo iff we define 



G 



ab 1 ...b 2 i\c 1 ...c N _2i- 



1 — V ) G |afe 1 



2\2N n {0) 

\ab 1 ...b 2 i\c 1 ...c N _2i- 1 \ 



(B.12) 



At - JV (27re M 2 ) 2iv 



Multiplication of ( P-4|) , with adapted indices and specification of Eg and A4, by 
thus leads to 



z(i+y)G aH 

Z- 1 G 



ab 1 ...b 2 i\c 1 ...c N _ 



A 



Gabt 



G 



Of/ \ 2 \^ n 1 1 \ P 6 i---M c i--- c iV-2i-i ( - r a6i...6 2 il c i--- c JV-2!-i G fc, 



|p[=Q 



+ 



64tt 2 \4(l + y) 
z 

G fel...fe2j-l|ci...C JV _2I-l 



A 



Gb 2j ...b 2l b 2 j-i G 



p — a 

b 2j ...b 2l a 



+ 



A 
A 

;i + ^) 2 



1 



JV-2J-1 

£ 

fe=l 



.62, 



^62^+1. .. 62(62j \ci...Cff-2l-l °"f>2j + l--- b 2ia|ci...Cjv_2!-l 

&2j - a 



C — c 

< - Jr c 1 ...c k b 1 ...b 2 ic k c k+1 ...c N _ 2 i_ 1 ^ r c 1 ...c k bi...b 2 iac k+1 ...c N _ 2l _ 1 

Cfc — a 



(B.13) 
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where p and \p\ are related by ( |5.2p . As always, for I = we have to replace G^ by G aa 
and G p6l by G pa . 

To solve the Carleman equation resulting from ( |B.13| ) for I > 1 we express w±y\ by 
( p.5| ) and use the same permutation trick as in section |6.2| . This fixes a and ci and reverses 
the order of b\ . . . b 2 i and c 2 . . . Cjv-2/-i, giving: 



G 



ab 1 ...b 2 i\c 1 ...c N _2i- 



A G{, 1 ...{, 2i _ 2a | Cl ... Cjv _ 2i _ 1 Gfc 2j ...fc 2; {, 2i _ 1 G6 1 ...;, 2j _ 2 b 2:| _ 1 | Cl ... Cjv _ 2i _ 1 Gfc 2j ...f )2ia 



^ + y?U (h-b 2l )(b 2j ^-a) 



A G;, 1 ...5 2j _ ia G; )2j . +liii ; )2; fe 2j | cl .. Xjv _ 2; _ 1 Gb 1 ...b2j-ib2jGb2j + i...b 2 ia\ci...c N ~2l-i 

Jl+Wj^ (h-b 2l )(b 2j -a) 

I ^ { ^ r c 1 ...c k _- l ab 1 ...b2ic k c k+1 ...c N _2i-i ^ c 1 ...c k _ 1 c k b 1 ...b 2 iac k+1 ...c N _2i- 1 , a 

+ (1 + ^) 2 ^ (6i-6a)(p*-a) ' 

In the simplest case (N — 4,1 — 1) the solution reads 

A G ai | c G a3Q2 G a2 | c G a3ai A G a2ai G Q3 | c G a2a3 G ai | c 



(1 + X> 2 (a 2 - a 3 )(a 2 - ai) (1 + X> 2 (a 2 - a 3 )(a 3 - ai) 
\ r 1 — c 

(1 + ^) 2 (a 2 -a 3 )(c-a 1 ) ' 1 ' J 

By general symmetry considerations, this solution should coincide with G c \ aia2a:j . A direct 
computation of G a \ clC2C3 would be similar to ( |B.11| ). 

B.4 Two cycles of even length: solution 

Our starting point is the equation ( |B.6| ) in which we rescale 

G Q 6 1 ...6 2i _ 1 | Cl ... Cjv _ 2! := ^"W^tf!..^^-* . (B.16) 

Accordingly, we multiply ( P-6|) by 11 Z {f+y)G b anc ^ obtain 

T ab 1 ...b2i-i\c 1 ...c N -2l 



Z ' G„ 



;i + y) G abl 

A 4 2 M 

= -Z(l + y)G abu .. b2l _ 1 { — (——-—) ^(|n|+l)G an , Cl ... Cjv _ 2! 

122:1=0 

A N ~ 2 



'ci...c k ac k ...c N _2i 



A / 4 \ 2 V^/i Ur P''i---f'2i-i|ci...c J v_2! < - r a6i...f) 2 ;-i|ci...c JV _2i G a6 



64tt 2 V^/i 2 (l + " ~ p-a 

I ^ \ ^ ■ '"J/ 1. ■■■"_■/ 1 '"J / 1 i ,\ it ■ 



A V V „ Ghr,-,* ...fo OI 1 /)o- 1^1 ...T.j\r of G, 



baj+i—6ai-i6aj|ci...cjv-a/ Ur &2j+i---f>2i-:Hci...c ] v_2i 

f>2j 
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i A Gb2 j+1 ...b 2 i- 1 b 2 j Gb 2j+1 ...b 2 i- 1 a 

(1 I W 2 / Cj fci-&2j|ei-CAr-ar 



1 I 3-') 9 ' ^ °1— 02j|Cl— CjV-2i ^ 



J 

N-21 



X) 



a 



, ^ t - J ci...c fc 6i...b 2 i-iCfeC fc+ i...c J v_2( ( - Jr ci...c fc 6i...6 2 i-iac fc+1 ...c JV _2i /o 1 ^ 

+ (^wh ~- ' ( ] 



where p and \p\ are related by ( |5.2| ). The thermodynamic limit 9 — > oo exists, and 
assuming Holder-continuity the usual steps (see ( |B.10|) ) lead to the following equation for 

G a b 1 ...b 2 i^ 1 \c 1 ...c N _ 2r 

^ a "T aG* g / V afe l- i '2i-l|ci...Civ_2i / ' rL a L LT »fel---fe2!-l|ci...C JV _2(J 

r A f A2 A ^ 2 I 

. A \^ rt Gb 2j+1 ...b 2 i^ 1 b 2j \c 1 ...c N ^ 2 i ~ Gb 2j+1 ...b 2l _ ia \ci...c N -2i 

\ t- 1 ri — C 

M I V\2 / / L 'b 1 ...b 2j \c 1 ...c N _ 2l - 



1 I 3^) 2 Ql...oaj|Ci...Cjy_2t ^ 



iV— 2£ 



A \ > G Cl ... Ck b 1 ...b 2l _ 1 c k c k+1 ...c N _ 2l G Cl ... Ck b 1 ...b 2l _ 1 ac k+1 ...c N _ 2 i / D 1t) \ 

+ <^y?h ~- ■ ( ' 

We first consider the case I > 2 where we use the permutation trick based on the invari- 
ance of Cr a 6 1 ...& 2l _ 1 | cl ..,c N _ 31 under simultaneous reversal of the order of bx . . .b 2 i-x and of 
C2 . . . cn-21- This eliminates the second line, and we obtain the algebraic solution 



G a b 1 ...b 2l ^ 1 \c 1 ...c N . 



21 



A ( , Gb 1 ...b 2j _ ia Gb 2j b 2j+1 ...b 2l _ 1 \c 1 ...c N _ 2l Gb 1 ...b 2j - 1 b 2j G a b 2j+1 ...b 2l _ 1 \ci...c N _ 2 i 

l+W^ (6sw-o)(6i-6a-i) 

Gbi...b 2j - 1 a\c 1 ...c N _ 2l Gb 2j b 2j+1 ...b 2l _ 1 ~ Gb 1 ...b 2j _ 1 b 2j \ci...c N _ 2 iGab 2 j +1 ...b 2l _^ 



(b 2j - a)(bx - b 2 i-x) 



\ N ~ 21 r - r 

A v { ^c 1 ...c k _ 1 ab 1 ...b 2 i_ 1 c k c k+1 ...c N _ 2l { ~Jc 1 ...c k _ 1 c k b 1 ...b 2l _ 1 ac k+1 ...c N _ 2 i / d1 m 

+ (ITf ^ ^ ' (&19) 



k=l 



This equation iteratively computes G 6 1 ...6 2 i_i|ci...c J v-2! fr° m the yet unknown function 
G a b\ Cl ...c N _ 2r Using the symmetry under exchange of the two cycles, the only unknown 
function is G a b\ clC2 which arises, for instance, in 



G 



abi&2&3|ciC2 



A Gb 1 aGb 2 b' i \cic 2 G b^b 2 G ab^\c\c 2 Gb 1 a\cic 2 Gb 2 bi Gfr^ b% \ c\ C2^a&3 

J+yy (6 2 _ a )( 6l _ 63 ) 

A I G ' nh\b 2 bzc\c 2 G Cl bib 2 bsac 2 G Cia bib 2 b^c 2 G ClC2 bib 2 b3a 



, " I auio 2 Q3Cic 2 ^ C1Q10203QC2 , ^ciaoiP203C2 ^ cic 2 o 1 D 2 osa \ /~r> \ 

+ (1 + ^) 2 V (a-aWx-h) + ( C2 _ a )( 6l _ 63 ) )■ 
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The yet undetermined function G a b\cd is the case / = 1 of ( |B.18| ): 
{- + ^- ) (aG ab]cd ) - —H a [ . G. b]cd \ 

= —Z(l + y)G a b(^r^I a \cd + jj-—yy{G ac( ic + Gadcdj^j 
\ a,, —a,, \ a — n „ 



A Gdcbc ~ Gdcba 



c — a 



7^ T—„ ' ( R21a ) 

(B.21b) 



(i + yy 
f A2 

Ia\cd '■= qdq G aq \cd ■ 

JO 

We regard I a \ c d as a function which is independent of G ab \ c d so that we can apply the 
Carleman formula ( |5.14|) . By linearity we first treat the auxiliary problem 



{- + ^ ) {aF abllci ... CN ) - — H a [. F. bllc ,.. CN \ 

A G C 2...CnC1 G C2 ^ CnC 



C\ — a 



(B.22) 



(1+^) 2 

By comparison with (|6.6| ) it follows that both sides of 

_v 

2 

G a b 1 c 1 ...c N — ^ ] Gb 1 c 1 ...C2i-i-^ab 1 \c 2 i- 1 ... Cjy (B.23) 

compute the functions F by iteration: 



satisfy the same Carleman equation. We can thus 

rp _ Gr a b lcl C2 

-Tafei|cic 2 — n i 

^6iciC2C3-^a&iC3C4 G a b 1 c\C2C'iCnG bia G bic\C2CzG ab\czc^ 

n, ~~ 7T, 



G 

rp , - T abiciC2C3C4 

F abi\ciC2C3C4 — n 

<J"6ici 



Gb x cxGb x 



C3 



(B.24) 



Since G a b is differentiable, these functions extend to coinciding indices q so that the 
inhomogeneity (1+ y )2 Gdcbc c ~^ dcba + (1+ A y)2 Gcdb <^ cdba in ( |B.21a| ) gives rise to the contribution 
F a b\cdcb + F ab \ dcdb in G ab \ cd . The in homo geneity proportional to Z(l + y)G ab needs to be 
treated by the Carleman formula Q5.14Q : 

G a b\cd = F ab \ cdc b + F ab \dcdb 

647r Ta^ (a) { Z(1 + y ^ Gah {-^ h \ cd + (1+ A y )2 ( G ^ + costf 6 (a) 

uAr,Q_l_ .a uA u. i 



+ e^Ui [e-*™ sin^(.)Z(l + ^ 



Fab\cdcb F ab \dcdb 

L$&(a) cosi9b(a 



A 

' M + (1 + y) 2 

A 



{G,cdc + C •dcd 



- ^ sin0&(a) cosi? 6 (a)GUZ(l + }>) (g^^ai + ^ ip^p (Gaafc+Gadcd)) 

- ^G ab Z(l + y)H K a [ Sln2 ^ ( ' } J. M + ^A^-^ + G^)) ] , (B.25) 
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where ( |5.19 ) has been used. We multiply by and integrate over b =: q to obtain an 
equation for I a \ cd := J^qdq G aq \ cd : 



X, 



a\cd 



1 



1 + y na J 
A 



qdq sin i? 9 (a) cos $ q (a)G a q \ + 7t 



X. 



\cd 



L 7T« 



q dq sin 2 $ 9 (•)£?, 



aq 



647T 2 



M 2 A 
J qdq {F aq \ cdcq + F aq \ dcdq ) + — 2 (C aC( j c +(j a rf C( j) 



X a \ cd :- Z(l + y) (-^Ia\cd + ^ + y j 2 ■ 
Using ^q = cot$ q (a) — cot$o( a ) an d inserting Z~ x from ( |5.7| ), we have 



A 



+ y) 2 

{Gacdc + G adcd )j 



(B.26a) 
(B.26b) 



r-X 



qdq sin$ 9 (a) cos $ q (a)G 



l + y na J 



aq 



A 



OqJ 



647T 2 



A ' G aq sm$ q (a) cos (# 9 (a) - # (a)) 



o 



All integrals in this equation, and also / q dq sin 2 $ q (b)G aq and / gd<7 (-F ag |cdcg + -Padded?) 

jo 
exist for A — >• oo, even in perturbation theory. Therefore, ( p.26a|) is again a Carleman 

singular integral equation (|5.13|) for X,\ cd where the input functions h and / only depend 

on the known data G ab and $&(a) and exist in the limit A — > oo. Its solution is given by 

( j5.14j ). Inserted into QB.25|) we thus get the (in principle explicit) solution 

647r A r647r 

G a b\cd = F ab \ cdcb + F ab \ dcdb - —— sin^ b (a) cos$ b (a)G ab X a \ cd - G ab H a — sin 2 $ b (»)X.\ cd . 

AO L a* 

(B.27) 

Together with ( |B.19| ) this gives the complete solution of all connected functions with two 
cycles of even length. 



sin$ (ci) 



(B.26c) 



B. 5 Remarks on functions with more than two cycles 

Although we are not going to work out the details, it is clear how to write down the 
Schwinger-Dyson equation for any planar S-cycle Appoint function and how to achieve 
its solution. By analogy with the (l+l)-point function and the (2+2)-point function we 
expect that the solution for the (Ni+ . . . +A£)-point functions with all iVj < 2 requires 
a case by case Carleman solution for known input data. If at least one Nj > 3, then 
the permutation trick gives an explicit algebraic formula for the (A^i+ . . . +Ag)-point 
function. 



C Perturbative expansion 

It is interesting to expand our solutions for the correlation functions into power series in A. 
These series should reproduce the expansion of the partition function into ribbon graphs, 
and indeed this agreement was for us a non-trivial consistency check of our equations. 
The Feynman rules are: 
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i + (i + y)(o+&) ^ or a (^) nne separating faces with indices a, b 

(—Z 2 X) for a (ribbon) vertex with four outgoing ribbons and index conservation at 



every corner, 

\2 



pdp 



(1 + y) / for every closed face with index p. 



o 



64vr 



These rules lead to a (for A — > oo) divergent one-particle irreducible (1PI) two-point func- 
tion r a b which therefore needs renormalisation through subtraction of the constant and 
linear Taylor terms. There is no subtraction for the four-point function: As result of the 
vanishing /3-function, the divergences in the ribbon graphs of the 1PI four-point function 
for A — )■ oo cancels exactly with the divergence of the wavefunction renormalisation Z. 

All of the following results are given for the limit A — > oo. Starting point is the master 
equation ( |5.38| ) for Gbo which yields G a o = ^-j- + 0(A). To next order we thus have 

G a0 = - ^ log(l + a) 

u 1 + a 64vr 2 (1 + a) K ' 

For the Hilbert transform we have 

A . w?[Grt] = __A_^W +oOT . 



64tt " 1 64tt 2 1 + a 
This yields for ( |5.16b| ) the expansions 

= A an _ /^A y o7r((l + a)log(l + a) - a log a) 
v; 647r 2 l + a + 6 V647T 2 / (1 + a + 6) 2 v ; ' v ' 

njooLa 1 a/ooro 1 _ A ~(l+o) log(l + 6) - a log a 2 

TCN-^o N - ^ + o(\ ) , 

which inserted into (|5.19|) give 

r 1 A (l + a)log(l + a) + (l + 6)log(l + 6) 2 

Ga6 " TT^Tb ~ 64^ (l + a + 6) 2 + ° {X } (C ' 2) 

for the 2-point function. This result is reproduced by renormalised one-loop ribbon graphs 
if we take the following expansion of y (see ([OOP ) and Z{\ + y) = e^W (see (jTgp) 
into account (strictly speaking, Z is only defined for finite A): 

^=^ + 0(A 2 ), Z(l+y) = l + 0(A). (C.3) 

It is iteresting to compare this with the perturbative solution of the true master 
equation (|5.26| ). This leads to T a = 1 + a + O(X) and then in next order 

T a = 1 + a + -Aj ((1 + a) log(l + a)+a log a) + 0(X 2 ) = A- 1 , (C.4) 
6471^ V / 647T tan Wo (a) 

in agreement with ( |C.1|) . We see this as good indication that in the limit A — > oo the 
master equations ( |5.26| ) and ( |5.38| ) are equivalent. 
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From ( |6.11| ) we then obtain the 4-point function G a bcd = G a bGbcGcdGdaT a bcd up to one 
loop, with the 1PI contribution 

A a - (1 + a) log(l + a) - c + (1 + c) log(l + c) 



abed 



(-A) 1 



64-7T 2 a — c 

A 6 - (1 + 6) log(l + 6) - d + (1 + d) log(l + d) 



+ 0(A 3 ) (C.5) 



64tt 2 b-d 
that agrees with the ribbon graph calculation. For the 6-point function, the result of 
( {3.12 ) can be arranged as 

Gabcdef = G ab G bc G cd G de G e fGf a (r abcd G ad T de f a +T bcde G be T e f ab +T cde fG c fTf abc 

^abedef ) > 

(-A) 3 (a-c)(l+e) log(l+e) + (c-e)(l+o) log(l+o) + (e-a)(l+c) log(l+c) 



r 



abedef 



64vr 2 



(a — c)(c — e)(e — a) 



(_ A )3 (6-d)(l+/)lo g (l+/) + (rf-/)(l+6)log(l+6) + (/-6)(l+rf)log(l+d) 

(b-d)(d -/)(/- 6) 



64vr 2 

+ 0(A 4 ), (C.6) 

also in agreement with the ribbon graph calculation. 
For the (l+l)-point function ( |B.11| ) we need 

A /clogc + clog2 aloga + alog2^ 



sin 



0(«)G.J 



64tt 2 V 1 + 2c 
A / (1 + c) log(l + c) 



1 + 2a 



+ 



log 2 



64tt 2 V(1 + 2c)(1 + 6 + c) (1 + 2c)(1 + 26) 



(l + 6 + c)(l + 26) 
Now ( |B.11| ) gives after straightforward but lengthy computation 



Ga\c G aa G cc (j? aacc G ac -\- r a | c ^ 
A 



2 r (l+c)log(l+c) _ (l+a)log(l+a) 



1+c+c 



1+a+a 



log(2) 



+ 0(A 3 ), (C.7) 



a|c 64vr 2 V (c-a) (l + a + a)(l + c + c 

where r a | c agrees with the Feynman graph computation 

r f°°Pdp , 

a|c J 64tt 2 (p + a + l)(p + p + l)(p + c + l) ^ 1 

For the (2+2)-point function (|B.27| ) we first have to provide some intermediate results. 
From ( B.24|) and QC.6 ) we have 



X a \ cd — G 2 cd \ G 2 c T ac ,j, 



G a bG cd G dc Gba (y£ abcdG adGbcT deba + ^ bcdc G bc G bc T cbab + G bc T abedebj 

'ZgLffqdq sin 2 ^(.)G a9 ] = (9(A 4 ), the 
gebraically with solution 

(1 - log(l + a))) + T adcd G 2 ad {\ - ^(1 - log(l + a))) 
+ 0(X A ) , (C.8a) 



" ab\cdcb 

Next we compute X a \ cd from (|B.26a| ). Since "H^ 
equations ( [B.26a| + [b.26cD are to C(A 3 ) purely a. 

A 2 



A 



64tt 2 

G add aac G ac J aad Jaadd d aacc 
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Ja 



pdp 



-A 3 ) 



;_A) 3 a-c + (l + c)(log(l+c) -log(l+o)) 



Jn 



64tt 2 (l+a+p) 2 (l+c+p) 64tt 2 (a - c) 2 

pdp (— A 3 ) 

o 64^ (1 + a + p) 2 (l + c + p) 2 
_ (-A) 3 2c - 2a + (2 + a + c) log(l + a) - (2 + a + c) log(l + c) 
" 64vr 2 (a - c) 3 

From (|C.8a ) we obtain after some rearrangements 



(C.8b) 



(C.8c) 



tin 



64?r 
(-A) 3 

647T 2 



sin 2 ^(•)X.| cd 



(1 + c) log(l + c) + al og(a) (1 + b) log(l + 6) - (1 + c) log(l + c) 



1 + a + c 



[l + a + c)(6-c) 



1 + d) log(l + d) + a log(a) (1 + 6) log(l + b) - (1 + d) log(l + d) 



(1 + a + d) 2 (l + a + d)(6-d) 

— G 2 cd G a b[j ccbG ca + JddbGda + dfebcc + J&Md) + (9(A 4 ) . 

These results and r afec(fcfe = J cca + J bbd give for ( |B.27| ) 

G a b\cd — G a bG c dGd c Gba abcdG a dG bc^ dcba + ^ abdcG ac G bdT cdba + ^ bcdcG b c G bc^ cbab 

+ FbdcdGbdGbdF dbab 4" r acdcG acG a<T caba dafea T J- a6|cd 



(C.8d) 



GadJaac + GbdJbbc + G ac J aa d + GbcJbbd + G c bJ C ca + G ca J cc b + GdbJdda + GdaJddb 



+ d aacc + daadd + J&& cc + Jf,Md + C*(A 4 ) . 



(C.9) 
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